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Excitation of 3 1 ' 3 /* and 3 1 ' 3 !) states of helium from the 
ground l 1 ^ state by electrons and positrons 



SURBHI VERMA and RAJESH SRTVASTAVA 

Department of Physics, University of Roorkee, Roorkee 247 667, India 

MS received 25 September 1997; revised 2 April 1998 

Abstract. A systematic theoretical study is carried out for the electron and positron impact 
excitation of the 3 ! ' 3 P and 3 li3 D states of helium from the ground 1*5 state. The results are reported 
for the differential cross sections, Stokes parameters and alignment and orientation parameters. To 
calculate the scattering amplitude which relate these parameters, a distorted wave approximation 
theory with different choices of distortion potential has been used. The results are compared with 
available theoretical and experimental data in the energy range of 40- 100 eV. 

Keywords. Excitation; electrons; positrons; helium. 
PACS No. 34.80 

1. Introduction 

Most of the pioneering experimental studies of the fundamental processes in electron- 
atom scattering such as elastic, excitation and ionization have been carried out with 
helium as target. Theoretically also any approximate method for a particular type of 
collision process before being recommended for the general use/application (i.e for 
heavier atoms) is first tested for hydrogen and helium. The studies on helium are more 
challenging due to the presence of two electrons in it which permits inter-electron 
repulsion to be introduced in a realistic manner. Further, helium being an inert gas is 
much easier to handle experimentally and because of this reason measurements with 
helium are more readily available than for hydrogen. In addition, the study of electron- 
helium scattering processes have wide ranging practical applications in the field of 
plasmas and lasers for their understanding. Also, helium is an abundant element and, 
therefore, helium emissions are widely observed in a variety of stellar and interstellar 
media, and are in turn used for determining abundances of this element. It is, therefore, 
important to provide quantitative electron collision cross sections for this target as well as 
to measure/calculate those scattering parameters which provide best test of a scattering 
theory. 

In this paper we focus our attention and deal with the electron (positron)-impact 
excitation of the n = 3 states of helium from its ground state with the point of view of 
differential cross sections and electron-photon coincidence results. Though electron 



ments as compared to those forn = 2 states which are order of magnitude higher. Also 
the close proximity of the n = 3 states make resolution of individual 3 3 ), 3 J D and 3 1 P 
states impossible even with the best available electron energy resolution. This can be 
understood from the fact that the 3 1 D and 3 3 D states are separated by less than 1 meV and 
the 3 1 P state is only separated from them by = 13meV. For example, differential 
scattering electron-energy-loss-spectroscopy (EELS) has yielded a significant quantity of 
data concerning both elastic and inelastic electron scattering channels. Using this 
conventional EELS, the measurement of differential cross sections for the excitation of all 
the strong 1*5 > n l P transitions in helium in principle should be possible, but have been 
limited to the energetically resolvable I 1 5 2 1 P transition providing for it the most 
accurate inelastic differential cross section (DCS) measurements. Until recently, (see for 
example paper of Khakoo et al [1]), no absolute experimental differential cross sections 
were available for the excitation of the n l P levels (n > 2) and the reasons being that 
conventional EELS from gaseous targets (employing electron beams of intensity of 
several nA) cannot provide the high energy resolution (< 13 meV) required to resolve the 
3 } P and 4 1 P levels from close-by levels, as well as providing adequate scattered electron 
signal for precise DCS measurements. 

However, in recent years, progress in the electron-photon coincidence angular and 
polarization techniques has greatly helped and provided the study of the excitation of the 
individual 3 1>3 P and 3 1>3 D states. As the scattered electron-decay photon coincidence 
technique takes the advantage of the fact that photons from the decay of these states, 
having different life times, can be readily isolated. Thus we find that most of the work 
published for individual n = 3 excitations in helium have become available only during 
the last couple of years when the coincidence experiments to study such excitations have 
started. To compare with these experiments many theoretical calculations such as R matrix 
(for low energies), close-coupling, first order many body theory and distorted wave Born 
approximation have been reported which specifically gave the results only for the range 
of scattering angles and the incident energies where the experimental data are available. 

In our previous papers [2, 3] we considered in detail the study of the n = 3 excitations 
in hydrogen and in continuation we give here our systematic similar detailed study of the 
n = 3 excitations in helium. We consider excitation of all the four 3 1)3 P and 3 1)3 Z) states 
from the ground I 1 5 state by impact of both the electrons and positrons. Note that with 
the positron impact only excitation of the 3 1 P and 3 1 D states will be possible. For the 
present study we extend our distorted wave approximation (DWA) calculation as used for 
hydrogen and helium [2-4]. Though the method used here is similar to the distorted wave 
Born approximation (DWBA) method used earlier by Madison and co-workers [6-8] and 
Katiyar and Srivastava [9], we have considered and tested now some extra features such 
as the incorporation of exchange and polarization effects in the distortion potential. Also 
as the different theories and experiments have reported results in a scattered manner for 
different parameters i.e. only at selected energies and (or) in particular range of scattering 
angles we present systematic detailed calculations for differential cross sections, Stokes 
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semi empirical attempt to obtain the DCS for the n = 3 levels using electron-energy-loss- 
spectroscopy and used first order many body theory (FOMBT) to calculate the 
contribution, to the signal, from the 1 1 S-3 1 D and 3 3 D. They gave results at incident 
electron energies of 29.7, 39.2 and 40 eV. Trajmar and co-workers [11] reported results at 
30, 50, lOOeV for the combined 3 1 P + 3 1 D + 3 3 D states alongwith giving an estimation 
for the DCS for the excitation of the 3 1 ? level at 80 and lOOeV. Allan [12] reported 
results in the low energy region i.e. 4eV above threshold. Further, Hummer and Burns 
[13] carried out the only direct relative differential cross section measurements for the 
excitation of the 3 1 P (with magnetic substate M/ = 0) state. To compare with these 
experimental data several theoretical calculations have been attempted [11,1416]. 
Madison [8] using the distorted wave Born approximation (DWBA) reported the DCS 
results at 80 and 100 eV. Cartwright et al [1 1] carried out the first order many body theory 
(FOMBT) calculations and reported the DCS for individual 3 1 P, 3 1 D and 3 3 Z> states at 
30, 40, 50 and lOOeVand thus also gave the DCS for combined 3 1 P + 3 1 D + 3 3 D states. 
In a separate paper Csanak et al [14] reported the DCS for magnetic substates with 
Mf = 0, 1 for the excitation of the 3 1 ? level at 50, 70, 80 and 100 eV. In this paper they 
also tabulated the experimental DCS data for different magnetic sublevels obtained using 
previously available data for DCS summed over different magnetic sublevels and the A 
parameter. Recently Fursa and Bray [15] carried out a convergent close coupling (CCC) 
calculation for the DCS for combined He(3 ! P + 3 1 D + 3 3 D) excitations in the energy 
region of 30-500 eV. In the low energy region Fon et al ([16], and references therein) 
carried out 19-state and 29-state jR-matrix calculations for the DCS of the excitation 
to individual 3 1 ' 3 / 3 and 3 1)3 D states of helium. In addition to these various reports of 
differential cross section results, the measurements and theoretical calculations have also 
been reported for the electron-photon coincidence studies on the excitation of the 3 1 P 
state of helium. Eminyan et al [17] obtained experimental data for the A and x parameters 
at incident electron energies between 50-1 50 eV whereas McAdams and Williams [18] 
reported the A and x parameters at an energy of 81.2 eV. In a later paper Beijers et al [1] 
measured the alignment and orientation parameters for the 3 1 P state of helium in the 
energy range 50-80 eV. Theoretical calculations for these parameters are those using 
FOMBT [19], DWBA [7] and CCC [15] methods. Further, in the low energy region Neill 
and Crowe [20] and Neill et al [21] reported the electron-photon coincidence parameters 
at 26.5 and 29.6 eV energies respectively. To compare with these low energy measure- 
ments Fon et al [22] carried out the /?-matrix calculations for A and x parameters. 

For the 3 3 P excitation of helium from the ground 1 1 S state some measurements for 
DCS have been carried out in the intermediate energy range viz. at 39.7 eV by Chutijian 
and Thomas [10], at 30, 50 and 100 eV by Trajmar et al [23] and at 40, 60 and 80 eV by 
Batelaan et al [24]. The results of Batelaan et al [24] were extracted from the existing 
data for electron-photon coincidence experiments. To compare with these available data 
in this intermediate energy range theoretical calculations using the FOMBT [23], DWBA 
[6] and CCC [15] have been performed. Further, in the low energy range Allan [12] 



Amongst these are the experimental data of Silim et al [25] and Beijers et al [26] at 60 eV 
for Stokes parameters and alignment and orientation parameters, Humphrey et al [27] for 
Stokes parameters at 39.7 and 23.5 eV and Donnelley et al [28] for the angular correlation 
studies at 29.6 eV. The theoretical calculations for these parameters have been carried out 
using FOMBT [29] and DWBA [6] in the intermediate energy region and using CCC [15] 
and /?-matrix [22] for low energies. 

For the He (l 1 ^ 3 1>3 D) excitations the only experimental data reported for DCS are 
due to Batelaan et al [24] who extracted the results from the available electron-photon 
coincidence experimental data as they did for 3 3 P excitation. The existing theoretical 
calculations for 3 3il Z) excitations in the intermediate energy range of 40-1 00 eV for the 
DCS are again from FOMBT [11], DWBA [6], CCC [30,31] and multi channel eikonal 
theory (MET) [32, 33]. As compared to the differential cross section measurements rela- 
tively more attention has been devoted to the experimental determination of electron- 
photon coincidence parameters for 3 li3 D excitations. Thus the work that started from the 
angular correlation studies of the 3 1 D states of helium (van Linden van den Heuvell et al 
[34]) has given impetus to many other studies and now various results for excitation of the 
3 1)3 D states in helium have been reported with increasing pace. The recent experimental 
papers [35-37] alongwith the reviews of Andersen et al [38] and Slevin and Chwirot [39] 
describe in a comprehensive manner the details of the angular and polarization correla- 
tion data available for these two transitions at different selected energies. Theoretically 
electron-photon coincidence results for the He (3 ! D) excitation process are obtained using 
the DWBA [6], FOMBT [40,41], MET [32, 33], CCC [30] and in the low energy region 
from /J-matrix ([16] and references therein). Similarly the He(3 3 Z>) excitation is also 
considered by Bartschat and Madison [6], Bray et al [31] and Fon et al [16]. 

2. Theoretical considerations 

2.1 T-Matrix 

The transition matrix for the electron (positron) impact excitation of helium from its 
initial state i to a final magnetic substate/M in the distorted wave approximation (DWA) 
can be written as 

T yu = <Xf u \V-U f \Axt), (1) 

where A is the antisymmetrization operator which takes into account in case of electron 
impact excitation the contribution of electron exchange between projectile and target 
electrons, V is the total interaction potential and is given by 



ri-r 3 



where ^/(/M) i s me initial (final) state wavefunction of the target helium atom and 
S,-(/)(l,2;3) is the initial (final) state spin function for the composite system consisting 
of the incident projectile and the target. F^i^ represents the initial (final) channel 
projectile distorted wave with the wave vector k,(k/) and the associated superscript +( ) 
indicates the usual outgoing (incoming) wave boundary condition. These distorted waves 
are the solution of 

[V* + k? w - 2t/,. (/) (r 3 )]F + H(k ( . (/) , rs) = . (4) 

Here U^ is the distortion potential in the initial (final) channel. 

On substituting the expressions for x}~ and Xf M fr m (3) into the r-matrix (eq. (1)) 
we get 

3)\V^ 



;l)}, (5) 

where <5z',-i is the Kronecker delta symbol and it goes to zero for positron impact causing 
the exchange term to be dropped from the T-matrix. 

Further, the scattering amplitude a M in the 'collision reference frame' is related to the 
DWA transition matrix element 7^ M for the electron (positron) impact excitation of 
helium from an initial state i to a final magnetic substate /M by 



The expression for the T-matrix (eq. (5)) can be simplified further for a specific transition 
by carrying out the integrations over the spin co-ordinates. In this paper we consider 
two types of transitions, namely the excitation from the ground 1 1 S state to excited 3 } D 
and 3 1 ? state i.e. singlet to singlet (SS) transition and the other is excitation from the 
ground 1 1 S state to 3 3 P and 3 3 D states i.e. singlet to triplet (ST) excitations. We give 
below a brief outline how to evaluate the r-matrix for each SS and ST excitation 
processes. 

2.2 Singlet-singlet excitations (l 1 ^ - 3 1 P and 3 1 D) 

In the singlet-singlet excitations the scattering takes place in the doublet mode i.e. the 
total spin of the system S = 1/2. The doublet mode spin function 5 1 ,- (S/) in terms of the 
usual Dirac matrices a and (3 for the composite system is given by [42] 



,2; 3)= -=a 3 (ai/3 2 - o&A). (7) 

V2 

Snhsritntintr ,^,/^\ in p.n (5\ anrl r.arrvincr out the; snin integrations minp the orthogonal 



<%, (8) 

<a,-|#>=0, 

we get 

T ifH = T* M -8 z ^T. (9) 

In equation (9) T^ u and 7|* are the spin averaged direct and exchange 7-matrices 
expressed as 

^M = / F "*( k /' r 3)^M( r i' r 2)(^-^(^))^ + (k,r3M-(ri,r 2 )dr 1 dr 2 dr 3) 

(10) 
T V = /^(k/,r 3 )^ M (r 1 ,r 2 )(V- t/ / C-3)) J F + (k J - ) r 1 )^(r 3) r 2 )dr 1 dr 2 dr 3 . 

(11) 
Thus the differential cross section for the singlet-singlet SS excitations is given by 



do 



_ \ A / rnnd c T ex 1 2-1 /n\ 

-4^Z^T.^ T ifM~^,-\T ifM \ \. (12) 

SS ^ M ' 



2.3 Singlet-triplet excitations (1 1 S - 3 3 P and 3 3 D) 

In the singlet-triplet excitations the scattering takes place in the doublet mode, i.e. in this 
mode also the total spin of the composite system S = 1/2. This is because of the 
conservation of the total spin of the system during the excitation process. The expressions 
for Si and Sf in the ST excitation are given by 

(13) 



S f =~ 7 = [2a ia2 fo - ^(a { j3 2 + a^}]. (14) 

vo 

On substituting these in the jT-matrix (equation (5)) and carrying out the spin integrations 
we get 



T ifM = x3<!>z',-ir-. ' (15) 

Here, the direct transition matrix 7^ M does not appear because the ST excitation cannot 
take place without an interchange of the projectile and target electron. Also note that ST 
excitations are therefore not possible by positron impact. 

Consequently the differential cross section for the ST excitation in helium atom is 
given by 

I da I 1 v-^ fy , /- _ ,o 



variational wave functions derived and used by Winter and Lin [43]. The ground state 1 1 S 
wave function for helium is given by 



(17) 
with 

Mo(r)=/?o(r)Koo(r), (18) 

where 

R (r)=Ae- a > r + Be- a * r . (19) 

The constant A, B, ot\ and a 2 can be obtained from Winter and Lin [43]. Further, the n 1>3 L 
states are having the following form 

0(n 1 ' 3 LM|n,r 2 ) = { M2 (n 1 ' 3 LM|r 1 )M 1 (Z|r 2 )+Mi(Z|r 1 )M 2 (n 1 ' 3 LM|r 2 )} 
V2 

(20) 
with 

m(Z|r) = J 1 (r)y 00 (r) (21) 

and 

u 2 (n l ' 2 LM\r) = /? 2 (n 1 ' 3 L|r)7 LM (r), (22) 

where 

^(r) = 2Z 3/2 exp(-Zr) (23) 

and 

n-L 

1 , (24) 



with Z = 2. The parameters b, c and N occurring in 2 have been taken from the paper of 
Winter and Lin [43]. The helium wavefunctions 0(n 3 L|ri,r2) for all n are orthonor- 
malized to each other and to the ground state 0(l 1 5'|ri,r2). Further all the one-electron 
wavefunctions M 2 (rc 1)3 L|r) have been set orthonormalized to one another and to the wave 
functions UQ(T) and wi(Z|r). 

2.5 Choice of distorting potential 

For the calculation of the distorted waves (using eq. (4)), the distortion potential U^ is 
taken as the sum of the spherical average of static potential (<j) a \V\(f) a } , exchange 
distortion potential and the polarization potential i.e. 

^(/ = y r + ^ I -i + yP 1 (with = * or/). (25) 



i-2 



(for I 1 S state) 



/ x M 

max(r 1 ,r 3 ) 



(26) 



-L 



o max(r 2 ,r 3 ) 



(for n 1 ' 3 ! state). 



is the exchange distortion potential in the initial (final) channel and is taken to be 
that of Furness and McCarthy [44] 

vgf fo)= H(i*?tn - V stat ( r3 - Ki*?</) ~ v * ta V3)) 2 - STrrp^)] 1 / 2 }. 

(27) 

Here p a (r) = f \(f> a (r, r')| 2 dr' is the spherical average of the one electron charge density of 
the helium atom in the initial or final channel. The value of the parameter r depends on 
the total spin of the colliding system and is taken from the paper of Vucic et al [45]. For 
the channel e He ( 3 L) it takes -1 for the doublet system and +2 for the quartet one. 
For e+ He (n { L), r = 1. To obtain V fl stat and p a we use the same target wave function as 
for the evaluation of the T-matrix. V?,l is the polarization potential and is taken as the 
widely used simple Buckingham type polarization potential proposed by Jhanwar and 
Khare [46] having the form 



(28) 



where a.& is the static dipole polarizability of the atom in the j'th state (see Miller and 
Bedersen [47]). The parameter d is taken from the paper of Mittleman and Watson [48] 
and is given by 

d = (-ad0 2 Z~ 1 / 3 } 1 ^ 4 (29) 

with /3 = I. Usually is considered as an adjustable (energy dependent) parameter which 
is determined by comparing with experimental data. Many workers (see e.g. Furness and 
McCarthy [44] and Verma et al [4]) took the adequate value of such that the experi- 
mental cross section for small scattering angles (where polarization is important) can be 
produced fairly accurately. Since the experimental data were not sufficiently available for 
all the excitations and to be consistent we took to be energy independent and equal to 
unity as recommended by Mittleman and Watson [48], 

Similar to our previous papers [3-5] three different choices of the distortion potential 
are made for the present study also, that is the distortion potential of either the ground 
state (II) or excited state (FF) is taken in both the initial and final channels, or the distor- 
tion nf thp. initial state, in the. initial r.hanne.l anH nf the. p.xnited state, in the final channel 



2.O blokes parameters ana alignment ana onentation parameters 

Blum [49] gave explicit expressions for the Stokes intensities in terms of the state 
multipoles (r(L) Q ) for the photons emitted in direction (0 7 , 7 ) in coincidence with the 
scattered electrons after impact excitation of an atomic state with orbital angular 
momentum L and an electric dipole transition to a state with orbital angular momentum 
L in the 'collision frame of reference' as 



1 1 2 
L L Li 



7 



3(2L+1) 1/2 



x G 2 (L){T(L)J 2 ) x sin 2 7 cos20 7 - G 2 (L)(7(L) 2 f ] )sin2^ 7 cos^ 
, G 2 (L)(r(L) 2 + )(3co S 2 7 -l)' 



(30) 



x [G 2 (L)(r(L) 2 f 2 )(l + cos 2 7 )cos2< 7 + G 2 (L)(r(L) 2 f ] ) sin20 7 cos0 7 
)sin 2 7 ], (31) 



x [G 2 (L){r(L) 2 f 2 }2cos0 7 sin20 7 +G 2 (L){r(L)+}2sin^ 7 sin0 7 ], (32). 



(33) 

where C(w) = (u//2-7rc 3 ) is a constant containing the frequency u of the transition, 7 is 
the decay constant and [(I^HrllL)! 2 is related to the oscillator strength of the radiative 
transition. GK(L) are the depolarization perturbation coefficients which account for the 
fine structure of the final excited state and are given by 



(34) 



where S is the electronic spin and / = L -f 5 is the total angular momentum of the atom. 
The G K (L] are normalized such that G (L) = 1 for all L. 
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/ 



M'M 

where 

<a(M')a(M)*) = * + (2S+l)a(M') s a(My*. (36) 

Here a(M} s is the scattering amplitude for excitation of the magnetic substate M in the 
channel with total spin S and 5,- is the spin of the initial target state. From the properties of 
the 3j symbol the values of rank K and component Q of the multipoles are restricted as 
follows: 

K < 2L and - Q < K < Q. 

The state multipole (T(L) J,} is a measure of the overall population of the atomic state and 
can be expressed in terms of the differential cross section a as 

>= ' (37) 



Further we see from equations (30-33) that (T(L}^ Q } describes the orientation of the state 
and thus the circular polarization of the emitted photons and the state multipole (T(L}^ Q ) 
describes the alignment of the atomic state and thus gives the anisotropic part of the 
coincidence rate and the linear polarization of the emitted photon. Also note that for 
dipole transitions it is impossible to achieve state multipoles with rank K > 2 since a 
single photon has an angular momentum of unity. 

The differential Stokes parameters for scattering of unpolarized atoms by unpolarized 
electrons P, (/ = 1-3) measured perpendicular to the scattering plane (0 7 = 7 = 90) are 

7(0) -7(90) 

1 7(0) +7(90)' ( j 



P -7(135) 

2 7(45) +7(135)' ( } 
7(RHC)-7(LHC) 

3 7(RHC)+7(LHC)' ( } 

Further, the Stokes parameter P$ is measured parallel to the scattering plane (0 7 = 90, 
</> 7 = 0) and is given by 

^7(0) -7(90) 

4 7(0) +7(90)' l j 

where 7(o: ) is the intensity of the light with the polarization vector in the a. direction with 
respect to the incident projectile direction and 7 (RHC) and 7 (LHC) are, respectively, the 
intensities of the right and left circularly polarized light components. 
Thus using equations (30-33) the Stokes parameters P,- (f = 1, 2 and 3) for any L > LI 



p 



1 1 2 
L L L, 



(42) 
(43) 



(44) 



and the Stokes parameter P^(9^ = 90, </> 7 = 0) for the same decay is expressed as 



(45) 



where 



P =- 



1 2 
L L, 



(46) 



and 



/*= 



1 1 2 
L L Li 

+ 



L 



(T(L}^}-G 2 (L}(T(L}+ 2 } . 



(47) 



Further, to describe the nascent charge cloud of the atom immediately after excitation 
Andersen et al [38] defined the parameters which can be related to the Stokes parameters. 
These are the alignment angle, 7, linear polarization, P/, height of charge cloud, poo, and 
the angular momentum transferred perpendicular to the scattering plane, Lj_. In this paper 
we have considered excitation from the ground l 1 ^ state to higher 3 1>3 P and 3 1)3 D states 
of helium which decay by photon emission to lower allowed 2 1)3 .S' and 2 1)3 P states 
respectively and for such cases the parameters defined by Andersen et al [38] are related 
to the Stokes parameters by the following relations 

(48) 
(49) 

(50) 



= -P 3 



2L 

2 + K 



and 



Here PI, P 2 , PS, and P^, are the reduced Stokes parameters which may be obtained by 
using the same relations as for measured Stokes parameters PI, PI, PS and P$ given by 
equations (42-47) but with all Gjt(L) taken as unity. 

For the sake of clarity we also give below the simplified explicit expressions of the 
Stokes parameters for the various excitations considered in this paper. 

2.6.1 For 3 li3 P excitation: The excited 3*P and 3 3 P states after excitation decay by 
dipole transition by photons to the lower S states viz. to 2 1 S and 2 3 S states. In this case 
L = 1 and Z/2 = 0, therefore, the Stokes parameters are given as 






(53) 

(54) 
(55) 



Here 



P = -^ G <r(3P) +> + i f-% <T(3P) 2 + ) + G 2 <r(3P)+}j (56) 

3V3 3 LV6 J 

and the depolarization coefficients are given by 

GO = GI = G 2 = 1 for 3 1 P state of helium, 

Go - 1, GI = 1, G 2 = 5! for 3 3 P state of helium. 

In addition to the above defined parameters (eqs (48-54)), Donnelly et al [28] and 
Humphrey et al [27] for 3 3 P excitation have also reported experimental data for the 
angular correlation function N i.e. the angular distribution of the coincidence photons 
which are collected in collision scattering plane. The expression for which can be 
obtained from equation (34) by taking </> 7 = 0. However, in terms of the well known A 
and x parameters after allowing for the depolarization of the decay radiation due to the 
fine-structure interaction among the triplet states this can be written as 



(57) 

where 

X = (al/(al + 2a\}} (58) 

and 



^t/j me oiuKcs parameier lor a ' u - 
process can be expressed in terms of state multipoles (T(3D)^ Q ) as 



ae-exciianon 



A/|G 2 (r(3D) 2 + )+G 2 <T(3D)J 2 
v < 
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where 



and 



such that 



I* = 



1U 



- G 2 (T(3D) 2 + 2 ) 



(60) 

(61) 
(62) 

(63) 
(64) 
(65) 



Go = GI = G 2 = 1 for 3 ! D state of helium, 

GO = 1, GI = ff, G 2 = j^5 for 3 3 D state of helium. 



3. Results and discussion 

Using the distorted wave approximation method as described in earlier section the results 
are obtained for differential cross sections, Stokes parameters and alignment and 
orientation parameters for the 3 1 ' 3 /* and 3 1 ' 3 /) excitation in helium from the ground 
1 1 S state. We have carried out calculations for both electron and positron impact 
excitations in the energy range 40 to 100 eV. The DWA calculations are carried out in all 
the four versions i.e. n, IF and FF and IIP. Although as mentioned earlier there are 
number of calculations available we have compared our various results only 
with the largely available FOMBT and the accurate convergent close coupling (CCC) 
calculations for sake of clarity and consistency. We have not made here comparisons of 
our results available with the distorted wave Born approximation calculations of Madison 
and co-workers [6-8] and Katiyar and Srivastava [9] as our results are found to be within 



Since me individual excitation or me ^f, $'L> ana yu states coma not DC resoivea 
earlier the combined measurements of the differential cross sections were reported. In 
figure 1 we therefore present the DCS for the combined He (3 1 / 3 + 3'D + 3 3 D) excitation 
at incident electron energies of 40, 50, 60 and lOOeV. These are compared with the 
available CCC results of Fursa and Bray [15] and the FOMBT results of Cartwright et al 
[11] at 40, 50 and 1 00 eV electron energies. They are also compared with measurements 
of Chutijian and Thomas [10] at 40 eV and of Trajmar [1 1] at 50 eV and 100 eV. At 40 eV 
our DWA results obtained using the n and IIP versions alongwith FOMBT and CCC 
calculations are in good agreement with the experimental results only upto scattering 
angle of 60. At 50 eVour DWA results in the FF model are in over all the best agreement 
with the experimental data though the other DWA calculations (i.e. in n, IF, IIP models), 
CCC and FOMBT also agree well with the experimental data below scattering angle of 
50. At 60 eV only our calculations are available which show similar features as at 50 eV. 
Finally, at 100 eV being the high energy as expected all the calculations agree closely 
almost upto scattering angle of 120 with themselves and with experiment. 

In figure 2 we present our DWA results as compared with other results for the 
individual 3 1 / 3 excitation at 40 and 80 eV. At these energies Khakoo et al [1] have reported 
recently the experimental DCS by direct measurement. In this figure, we included the 
experimental data of Khakoo et al [1] and the CCC and FOMBT calculations. From the 
figure we see that at 40 eV our n and HP results are the closest to the experimental data 
while at 80 eV except our IF model all other models give results which are in reasonable 
agreement with available experimental and other theoretical results over the entire 
scattering angle range. 

In figure 3 the DCS results for excitation of the Mf = substate of the 3 1 P state are 
shown at 70 eV energy. These are compared with the directly measured relative DCS of 
Hummer and Burns [13] and the FOMBT calculations of Csanak et al [14]. Our DWA 
results obtained using the FF model are in reasonably good agreement with the experi- 
mental data. Note that the relative DCS of Hummer and Burns [13] have been norma- 
lized, in a manner similar to the one used by Csanak et al [14] to our DWA results in the 
FF model at scattering angle of 10. In figures 4, 5 the DCS results for the excitation of 
the different magnetic substates (i.e. Mf = 0, 1) of the 3 1 P state are presented at 50, 80 
and 1 00 eV energies. These are compared with the only available FOMBT calculations of 
Csanak et al [14] as well as with the estimated experimental data obtained utilizing the A 
parameters and the DCS as explained and tabulated in their paper. Comparing the results 
for Mf = in figure 4 we find our FF model results are in the best agreement with the 
experimental data at 50 eV whereas the FOMBT results lie in between our n and IIP 
calculations and differ from FF results. At 80 and lOOeV the different theoretical results 
have come closer and the experimental data lie in between FF and other calculations. In 
figure 5, the results are presented for DCS for 3 l P(Mf = 1) excitation. From this figure 
we see that our different DWA results at 50 eV compare well with the available experi- 
mental and FOMBT results below scattering angle of 60. Beyond 60 the experimental 
data tend to lie above the various theoretical results. Further, on comparing the different 
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Figure 2. Differential cross section in units of cm 2 sr~ 1 for the electron impact 1 1 S- 

3 1 P excitation in helium at , = (a) 40 eV and (b) 80 eV. present FF results; 

present IF results; present n results; present IIP results; 

-- FOMBT results of Cartwright et al [11]; --- CCC results of Fursa and Bray 
[15]; A experimental data of Khakoo et al [1]. 



results for Mf = and for Mf = 1 shows that their nature are entirely different which on 
combining finally show the behaviour of Mf = (see figures 1 and 4). 

In figures 6, 7 and 8 we present the Stokes parameters PI, P 2 and P 3 for the decay of 
photons from 3 } P state to 2 1 S state which we calculated using equations (53-56). In these 
figures the parameters of Andersen et al [38] (viz. P/, Lj_ and 7) for the excited 3 1 P state 
obtained using equations (48-50) are also shown. The results are reported at incident 
electron energy of 50, 80 and lOOeV as most of the experiments and other theoretical 
calculations have been carried out at these energies. 

In figure 6 we compare our various DWA results for Stokes parameters at 50 eV with 
the measurements of Eminyan et al [17]. The agreement between the experiment and 
theory is good for all Stokes parameters i.e. PI, PI and P 3 . As seen for the DCS results 
out of our different DWA results the FF model calculations are in the overall best 
agreement with experimental data. Further, the inclusion of the polarization potential in 
the distortion potential (HP) also improves the agreement with the experimental data. In 
this figure the displayed alignment and orientation parameters (viz. P/, L_ and 7) are 
compared with the CCC calculations [15] and the experimental data of Beijers et al [7] 
and Eminyan et al [17]. On comparison for Lj. results we find our DWA calculations 
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Figure 3. Differential cross section in units of cm 2 sr~ 1 

1 1 S-3 1 P(M/ = 0) excitation in helium at E t = 70 eV. 

present IF results; present II results; present IIP results; 

--- FOMBT results of Csanak et al [14]; experimental data of Hummer and 
Burns [13]. 



for the electron impact 
present FF results; 



the 7 parameter the CCC results alongwith the FF model are in better agreement with the 
experimental data as compared to other results. For the PI parameter, again the FF model 
can be said to work the best in terms of overall agreement with the experimental data 
alongwith CCC calculations, although the second minimum of the CCC calculation is 
absent in our DWA results. 

In figure 7 the Stokes parameters PI, Pa and PS at 80 eV are compared with the 
measurements of Eminyan et al [17] and Me Adams and Williams [18]. The agreement 
of the DWA results with the measurements of Eminyan et al [17] is good for all the PI, 
PI and P 3 Stokes parameters. Also the experimental data in the larger scattering region 
due to Me Adams and Williams [18] though agree well with our DWA calculations 
for 7*2 but lie below and above our results respectively for PI and P^, Stokes 
parameters. Further in this figure the results for Lj_ 5 7 and P; parameters are compared 
with the experimental data of Beijers et al [7], Eminyan et al [17] and McAdams and 
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Figure 4. Differential cross section in units of cm 2 sr" 1 for the electron impact 1 1 S- 

3 l P(M f = 0) excitation in helium at / = (a) 50 eV, (b) 80 eV, (c) lOOeV. 
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present EDP results; --- FOMBT results of Csanak et al [14]; estimated 
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Figure 6. Stokes parameters PI , PI and P 3 and alignment and orientation parameters 
PI, 7 and L for electron impact I 1 S-3 1 P excitation in helium at ,- = 50 eV. 

present FF results; present IF results; present EL 

results; present IIP results; --- CCC results of Fursa and Bray [15]; D 

experimental data of Eminyan et al [17]; experimental data of Beijers et al [7]. 



[19]. Features for these parameters are similar to figure 6. The FOMBT results are 
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Figure 7. Same as figure 6 but for E- t = 80 eV and ---, FOMBT results of 
Csanak and Cartwright [19]; A, experimental data of McAdams and Williams [18]. 



Through figures 9-12 results are presented for DCS, Stokes parameters and alignment 
and orientation parameters for the positron impact I 1 5 3 1 / 5 excitation in helium. These 
results are presented at positron energies of 50 and lOOeV. Since there are no 
experimental results to compare we have presented only our different results. 

Figure 9 shows the DCS in the different H, IF, FF and IIP models at 50 and lOOeV 
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Figure 8. Same as figure 6 but for E t = 100 eV and , FOMBT results of Csanak 

and Cartwright [19]. 



corresponding electron-impact results (see figure 1) we find that the results for positron- 
impact fall off quickly for larger angles as compared to the results for electron-impact 
excitation. This may be due to absence of exchange contribution in positron-impact 
results. Rest other features are similar to as for electron-impact case and the positron- 
impact results are comparable with the electron-impact results. From figure 10 for the 
3 ! P(M/ = 0) excitation, the behaviour is similar to the DCS summed over all magnetic 
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Figure 9. Differential cross section in units of cm 2 sr~ 1 for the positron impact 

1 1 S - 3 J P excitation in helium at ,- = (a) 50eVand (b) lOOeV. present 

FF results; present IF results; present II results; present 

IIP results. 



with the electron-impact results, we find that the minimum in the electron impact results 
is totally missing here (see figure 4). Thus the presence of the minimum may be said to 
occur because of the difference in the nature of interaction potentials in both the cases of 
electron and positron-impact. 

In figure 12 for the sake of illustration in case of positron impact, the results at positron 
impact energy of 50 eV are presented for Stokes parameters (Pi, PI and PS) and 
alignment and orientation parameters (Lj_, P; and 7). Again no other results are available 
for these parameters to compare with our DWA results. The inter-comparison of these 
results with the corresponding electron impact results (see figure 6) gives the most 
striking feature that in contrast to the electron impact results, the Pa, L and 7 parameters 
retain the same sign throughout the angular region. Also the sharp maxima and minima 
observed in the electron impact case are missing in the different results for positron 
impact excitation. 




30 60 90 120 150 

Scattering angle (deg) 



180 



10 



10"' 



10 



10 '" 



10 



30 60 90 120 150 
Scattering angle (deg) 



180 



Figure 10. Differential cross section in units of cm 2 sr" 1 for the positron impact 
1 1 S-3 1 P(M/ = 0) excitation in helium at E,- = (a) 50 eV and (b) lOOeV. Notation 
same as in figure 9. 
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Figure 12. Stokes parameters PI, PI and P 3 and alignment and orientation 
parameters P lt 7 and L L for positron impact 1 1 S-3 } P excitation in helium at 
E; = 50 eV. Notation same as in figure 10. 



nven. These are compared with the CCC [15] and FOMBT [23] calculations and with 
Experimental data of Chutijian and Thomas [10] Trajmar et al [23] and Batelaan et al 
. Since Batelaan et al [24] reported relative DCS, we normalized these as best fits to 
DWA results in the FF model. On comparing the results at 40 eV, the DWA and 
dBT results tend to lie higher than the experimental data. Also, the CCC results are 
tivelv in better agreement with the measurements in the entire angular range. Since 
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Figure 13. Differential cross section in units of cm 2 sr 1 for the electron impact 
1'5-3 3 P excitation in helium at E { = (a) 40 eV, (b) 50 eV, (c) 80 eV and (d) 100 eV. 

present FF results; present IF results; present n 

results; present HP results; -- FOMBT results of Trajmar et al [23]; --- 

CCC results of Fursa and Bray [15]; Q experimental data of Chutijian and Thomas 
[10]; experimental data of Trajmar et al [23]; A relative experimental data of 
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Figure 14. Differential cross section in units of cm 2 sr 1 for the electron impact 
1 1 S - 3 3 P (M f = 0) excitation in helium at E-, = (a) 50 eV, (b) 80 eV and (c) 100 eV. 
present FF results; - - present IF results; present II 
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Figure 15. Same as figure 14 but for electron impact 1 1 S-3 3 P (Mr = 1) excitation 



lent of our results with both the available experimental data alongwith theoretical 
itions is reasonably good. Next in figures 14 and 15 we present the magnetic 
el DCS for the 3 3 P excitation in helium at 50, 80 and 100 eV energies. Since no 
previous DCS data are available for the individual magnetic sublevel we display 
ur DWA results in all the four models (viz. II, IF, FF and IBP). For the 3 3 P(M/ = 0) 
ion the cross sections in the IF model only show a pronounced minimum as 
red with others. Further, the HP results lie lower showing that the inclusion of 
nation potential tends to lower the magnitude of DCS. An interesting feature that 
noticed for 3 3 P(M/ =1) excitation is the very small magnitude of DCS for small 
and then the appearance of maxima in the scattering angle range of 20-30 which 
ill off rapidly with the increase of scattering angle. Further, the maximum is shifted 
Is lower angles and the width of the maximum also decreases with the increase in 
'. Also note that in contrast to the singlet-singlet scattering where all the four 
s give nearly the same values of DCS at low scattering angles, for the singlet-triplet 
ing the results obtained by different models (viz. n, IF, FF and IIP) differ 
iably. The difference in the DCS results for the singlet-singlet (1 1 S 3 1 P) and 
t-triplet 1 1 S 3 3 P transitions can be due to reason that the singlet-triplet transition 
present system (e + He) is caused only by the electron exchange interaction, 
igure 16 we show the coincidence count rate (AT) calculated using equation (57) 
in figures 17 and 18 our DWA calculations respectively at incident electron energies 
1 60 eV are shown as compared with other available results for Stokes parameters 
>, PS and alignment and orientation parameters (L, 7 and P/). The results in 
16 show the coincidence count rate N as a function of photon angle 7 . We have 
. our DWA results compared at an energy of 40 eV and scattering angle 60 with the 
rements of Donnelly et at [28] and Humphrey et al [27] which have been carried 
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Figure 16. Coincidence count rate W as a function of photon angle for the electron 
impact 1 1 S-3 3 P excitation in helium at j = 40eV. Same as figure 14 and A 
experimental data of Donnelly et al [28]; Q experimental data of Humphrey etal [27]. 
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Figure 17. Stokes parameters P\, PI and P-$ and alignment and orientation 
parameters PI, 7 and Lj. for electron impact 1 1 S-3 3 P excitation in helium at 

EJ = 40 eV. Notation same as in figure 16 and FOMBT results of Csanak and 

Cartwright [29]. 



out at 39.7 eV energy and 60 scattering angle. From the figure we find that the different 
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Figure 18. Same as figure 3.18 but for J?,- = 60eV. Notation same as in figure 17 and 
--- FOMBT results of Csanak and Cartwright [29]; experimental data of Beijers 
et al [26]; A experimental data of Silim et al [25]. 

on energy of 39.7 eV. Our various DWA results qualitatively agree with the available 
imental data for PI and P 3 parameters although for PI the measurements of 
?hrey et al [27] at 90 are not in so good agreement with our results. Further for L, 
1 PI parameters our results are compared with the reported experimental data of 
Dhrey et al [27] and the FOMBT [29] results for Lj_ and 7. The agreement for LI 
ien the different theories and the experimental data is good whereas for the 
cnent angle 7, the experimental data above 40 lie much below all theoretical results 
i for Pi the depth of minima in the experimental results is not so sharp and is less 
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parameters are presented at 60 eV alongwith the data of Beijers et al [26] and Silim et al 
[25] as well as the FOMBT [29] results available for Lj_ and 7. Our DWA results agree 
well with the experimental data for all these parameters alongwith the FOMBT results. 
Further as expected with the increase in energy the various DWA results tend to come 
closer. 

3.3 3*D excitation 

In figure 19 the DCS results are presented for the 3 1 D excitation in helium at electron 
energies of 40, 60 and 80 eV. These are compared with the only available CCC 
calculations of Bray et al [30] and the FOMBT calculations of Cartwright et al [11] at 
40 eV. The relative DCS data of Batelaan et al [24] at 60 eV are also reported which we 
compared after normalizing them as best fit to our results of FF model. We see from the 
figure at 40 eV that all the theoretical results are in reasonable agreement among 
themselves except CCC results which are much larger at smaller angles. At 60 eV our 
DWA results are in agreement with the experiments of Batelaan et al [24] while rest of 
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Figure 20. Stokes parameters P\, P 2 , P$ and P 4 for the electron impact 

excitation in helium at E t = 40 eV. present FF results; present 

IF results; present n results; present IIP results; --- CCC results of 
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Figure 21. The alignment and orientation parameters PI, 7, Lj_ and poo for 
electron impact 1 1 S~3 1 D excitation in helium at EI = 40 eV. Notation is same a 
figure 20. 
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Figure 23. Same as figure 21 but for , = 60 eV. 
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Figure 25. Same as figure 21 but for ,- = 80 eV. 
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Figure 26. Differential cross section in units of cm 2 sr l for the positron impact 
1 1 S-3 1 D excitation in helium at E { = (a) 40 eV and (b) 80 eV. present FF 



helium. These are compared with the experimental data of Donnelly et al [36], Mikosza 
et al [37] and Batelaan et al [35] and the CCC results of Bray et al [30]. Our DWA results 
are in reasonably good agreement with the experimental data though for the P/ and 7 
parameters the experiments tend to lie qualitatively higher than our results. Further, as 
expected from other previous calculations using the distorted wave approximation our FF 
model results are found to give best agreement with the experimental data. Similar to the 
results presented in figures 20 and 21, we give in figures 22 and 23 our corresponding 
various DWA results at 60 eV electron energy which are compared with the only available 
experiments of Batelaan et al [35] and Mikosza et al [37]. The features observed are as at 
40 eV. Finally, in figures 24 and 25, similar results at SOeVare given and compared with 
the single available experimental data of Mikosza et al [37]. As expected with the 
increase in energy the different DWA results come closer and are in good agreement with 
experiment. As can be seen from all the figures 19-25 and comparing our calculations in 
different models we find that the results are quite sensitive to the choice of the distortion 
potentials in the initial and final channels. In addition, on comparing the results in the IIP 
and II models (i.e. with and without the polarization potential) it is noted that the 
inclusion of the polarization effect in the distortion potential improves the agreement with 
the experimental data slightly. 

In figure 26, the DCS results are presented for positron-impact 3 { D excitation at 40 and 
80 eV incident positron energy. Here the FF results are higher than the II and IF results 
whereas as seen for the electron-impact excitation the results in FF model were lowest 
with respect to the other models. 

3.4 3 3 D excitation 

In figure 27 the DCS results for 3 3 > excitation at 40, 60 and 80 eV are shown. Our DWA 
calculations in all the models are given. At 40 eV they are compared with the available 
CCC [31] and FOMBT [11] results. The relative DCS normalized as best fit to our FF 
model results of Batelaan [24] are also presented at this energy. We find the CCC results lie 
lower than the DWA results in the smaller angular region and fall rapidly in large angular 
region. The FOMBT results are similar to the DWA results in small scattering angles. At 60 
and 80 eV electron energies since no other results are available for comparison only our 
DWA results are presented which behave like the results at 40 eV. In figures 28 and 29, 
respectively the results for electron-impact excitation of the 3 3 D state are presented for the 
Stokes parameters (Pj, PI, P-$ and P4) and the parameters (7, P;, Lj_ and poo) m & U the 
different models of DWA. These are compared with the experimental results of Batelaan 
et al [35] and Crowe et al [31] and the CCC results as given in Crowe et al [31]. For this 
transition the CCC results are found to have better overall agreement with the experimental 
data than our DWA results. In figures 30 and 31 similar results for these parameters at 
60 eV are shown. Due to non-availability of any other measurement or calculation we have 
presented only our DWA results. The observed features are similar to those for 40 eV. 
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Figure 27. Differential cross section in units of cm 2 sr" 1 for the electron impact 

1 1 S-3 3 D excitation in helium at E t = (a) 40 eV, (b) 60eVand (c) 80 eV. 
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Figure 28. Stokes parameters PI, PI, Pj and P^ for the electron impact 1 1 S-3 3 D 

excitation in helium at / = 40 eV. present FF results; present 

IF results; present II results; present IIP results; --- CCC results of 

Bray et al [31]; D experimental data of Batelaan et al [35]; A experimental data of 
Crowe et al [31]. - 
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Figure 30. Stokes parameters PI, P 2 , P 3 and P 4 for the electron impact 1 1 S-3 3 D 

excitation in helium at ,- = 60 eV. present FF results; present 

IF results; present n results; present IIP results. 
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4. Conclusion 

An extensive theoretical study in the distorted-wave approximation has been carried out 
for the electron impact 3 li3 P, 3 ! ' 3 D excitations from the ground I 1 S state in helium. 
Results are reported for the Stokes parameters and alignment and orientation parameters 
and DCS in the energy range of 40-1 00 eV. These are compared with recent theoretical 
and experimental results. The DWA results are found sensitive to the choice of the distor- 
tion potential used in the calculations of the distorted waves. As well known through the 
comparison of first Born results and distorted wave results that the effect of distortion 
potential which is dominated by static and exchange potentials is reflected in the results 
for the differential cross sections and different parameters more at the intermediate and 
backward scattering angles except that the part of polarization potential effects the 
forward scattering only which is not very dominant. Further the results for Stokes 
parameters and alignment and orientation parameters are more sensitive in comparison to 
the DCS results to the model employed for the calculation and hence these provide a 
better test for the usefulness of the various models. The calculation has also been 
extended to positron-impact 1 1 S 3 J P, 3 ] D excitations. In contrast to the electron impact 
case, we see the striking feature, that PI, P^, PS and Lj_ retain the same sign throughout 
the entire angular range for positron-impact scattering. 
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Abstract. The dissociation energy of twelve diatomic molecules has been determined by fitting 
four-parameter potential function U(r) = D c [[l exp{ b(r r e )}]f[\ Cexp{ b(r r e )}}] 2 to the 
true Rydberg-Klein-Rees (RKR) curves for their fifteen electronic states using the mean square 
deviation as the criterion for the selection of the best fit. Average deviation of D e has been found to 
be 2.7% as compared to 20.5% obtained with Lippincott's potential function for these molecules. In 
addition the anharmonocity constant u e x e has also been calculated for the same electronic states 
yielding average mean deviation 8.9%. 

Keywords. Dissociation energy; four-parameter potential; anharmonicity constant. 
PACSNo. 33.15 

1. Introduction 

The spectroscopic data for most of the diatomic molecules is available for vibrational 
levels which are considerably below the dissociation limit. So different extrapolation 
methods have been employed to find the dissociation energy of these molecules and 
Birge-Sponer [1] approach is most commonly used for this purpose. Obviously, the results 
inherently suffer from an uncertainty determined by the extrapolation gap and the method 
of extrapolation. In addition, it has long been known [2] that vibrational energies near the 
dissociation limit and hence the vibrational spacings depend on the functional form of the 
long range potential. This, in turn, leads to a marked deviation of graphical extrapolation 
from the accurate estimate of the dissociation limit. Using asymptotic approximation for 
long range potentials of the form 



where D e is the dissociation limit of the potential and C n is given by 

T.pRnv nnH rn- workers r9-^ll> have develoned WKB based extraoolation techniaue to 



state 
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be pointed out that a proper RKR analysis yields a much more accurate estimate of 
D e [5]. 

In 1990, Wei Hua introduced a four-parameter (FP) potential function [6], 



U(r}=D e 



r e }} 



-fc(r - r e }} 



(3) 



where b = (I - C}a, and C is an adjustable parameter having \C\ < 1 with a= ^Jk/J2D e . 
The FP potential fits very closely to the RKR curves compared to Morse and other 
analytical functions [6], when the potential domain extends to near dissociation limit. It is, 
therefore, natural to expect that FP potential would predict the dissociation limit precisely. 
We have made an attempt to investigate this aspect of the FP potential and the results 
constitute the content of this communication. Wei Hua considered fifteen electronic states 
of twelve molecules to examine his FP potential and the same electronic states have been 
used in the present studies as the spectroscopic constants, the RKR potential curves and the 
dissociation energies are known very accurately for these states. The molecular constants 
and respective dissociation energies are given in table 1. 



2. Estimation of the dissociation energy 

Because of the anharmonocity the molecules in the highly excited vibrational levels spend 
most of their time in the neighborhood of outer turning points. Therefore, the vibrational 



..Jin--- 1 (4) 

where j stands for various potential functions, has minimum value. With the rough 
estimated values of D e from the linear extrapolation method, i.e., D e = u%/4(jj e x e and C 
from the relation [9] 



1 \ 120 37 2 2 17 G 
~ ~ 



where 

x = , F = -?-- , and G = 



for the whole potential curve as the starting points, the potential energies corresponding 
to FP potential have been obtained by changing D e and C in steps of 100 cm" 1 and 
0.1, respectively. The resulting values are correlated with the experimental results and 
the values of D e and C corresponding to minimum value of Yj are chosen for the second 
scan. In the second and subsequent scans, the values of D e and C are varied in steps of 
decreasing magnitudes. The process is continued till Yj is minimum. The last scan has 
been carried out .by changing D e by 1 cm" 1 and C by 0.001. The value of D e that 
corresponds to minimum deviation Yj in the final scan is regarded as the best estimate of 
D P and relevant C value. 



3. Results and discussion 

Dissociation energies of fifteen electronic states of molecules LJ2, Na2, K2, Rb2, Cs2, Cl2, 
12, IC1, HF, H2, CO, XeO have been obtained employing the molecular constants listed in 
table 1 and the results are presented in table 2. For comparison, values of D e obtained from 
Lippincott's analytical potential function [10] 



U(r) = D e (l - e-*){l - abx l ^QXp[-bx l/2 ]}, (6) 

where 



F = -TT- , a = ; rr 5 W = and 

65*' l + (5F/4)' r e \ab) 

x= (r e , b 1.065; 

which is considered to be best among the three parameter potential functions for the 
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have been computed for these electronic states using the relation 

u) e x e = r 2 (l + C + C 2 } = B e A(l + C+C 2 ) (7) 



obtained in our earlier work [9]. Here, we have used C calculated above for the outer 
limb. These results together with the experimental values of u> e x e and the ones found by 
us in ref. [9], are compiled in table 3. 

Earlier works on the studies of electronic states of the molecules subjected to the 
present analysis (see the references in table 1) show that the properly executed NDE 
method of LeRoy et al [2-4] and RKR approach yield accurate values of D e . So we have 
taken these D e values as experimental (being the best prediction) to compare our findings 
and have given these in the second column of the table 2. It may be mentioned that D e for 
HF included in this table has been obtained earlier from the limiting curve of dissociation. 

A comparison of our results constituting third column of table 2 with experimental 
values listed in column 2 shows that the extent of agreement depends on the method. For 
the electronic states Li 2 (X 1 S+), Na 2 (X 1 +), C1 2 (X 1 S+), I 2 (X 1 S+), ICICLE*) and 
Cs 2 (X 1 E+) for which experimental D e values have been found with NDE technique (that 
includes the long range forces) the difference is less than 2%. However, in A and A' 
electronic states of IC1 the deviations are respectively 2.5% and 2.6%. The deviation in 
the case of the electronic state A can be viewed in the light of reported vibrational 
perturbations near dissociation and the merging approach (merging the constants from 
individual absorption band and emission band through least square fit) used to obtain D e 
[16]. Similar explanation is valid for state A' also since there is a possibility of pertur- 



literature, it is found that the highest observed vibrational levels in K 2 , CO and Rb 2 viz. 
v = 37 at 63334.333 cm- 1 for CO [13,14], v = 50 at 3687.749 cm- 1 for K 2 [11], v = 72 
at 3313.143 cm" 1 for Rb 2 [12], giving E max /D e ~ 0.70, 0.83, and 0.84 respectively, where 
max denotes the energy of highest observed level, are fairly below the respective 
dissociation limits. That means, RKR potential curve has been drawn accurately only up 
to the observed levels. Moreover, since RKR does not provide any direct means of extra- 
polation beyond the observed levels [2], indirect methods have been employed to obtain 
D e values which are likely to be at the variance from the accurate values. In the case of 
XeO and H 2 the levels are observed only up to v = 8 and v = 14, respectively, so that. D e 
are expected to be more erroneous. 

Long range forces have dominant influence [2] on the vibrational levels close to 
dissociation limit (~ top 10% of the potential well). With a view to examine the perfor- 
mance of the FP potential function close to the dissociation limit, D e (denoted as D e 
(FPj/s) in table 2) and C (referred to as C m i n (l/3) in table 2) have been derived consi- 
dering only the highest one-third of the known vibrational levels. But for the CO 
electronic state, the results (given in column 10 of table 2) for all other states have been 
improved significantly. The slight adverse result of CO is understandable since the known 
data constitute only 70% of the potential well starting from equilibrium position. This 
observation, besides providing better accuracy of D e , also strengthens the earlier findings 
[6], that FP potential describes experimental potential curves more accurately especially 
when the potential domain extends to near the dissociation limit. 

Values of vibrational constant uj e x e , when computed from the outer limb of the RKR 
curve show a marked improvement over the ones calculated using inner wall as well as 
the outer wall (i.e. r min < r < r max ) of the potential curve (table 3). The average mean 
deviation from the experimental values has been reduced to 8.9% from 17% [9]. Further- 
more, since uj e x e has been calculated using outer limb of the potential curve, the findings 
are in conformity with the fact that, because of anharmonocity, the vibrational spacings 
and other properties of diatomic molecules [2, 7, 8] are explained mainly by the outer 
limb of the potential curve. 

4. Conclusion 

FP potential function of Wei Hua can be used to predict dissociation energy of the 
electronic states of diatomic molecules. The accuracy of prediction is close to the one 
obtained by NDE long range expansion technique or by RKR potential curve constructed 
up to near dissociation limit. Besides, FP potential function, through equation (7), predicts 
ijj e x e value better than the known three parameter analytical functions. 
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Abstract. Fluorescence spectrum of Dy 3+ , Dy 3+ -Ho 3+ , Tb 3+ and Tb^-Er 3 " 1 " doped in calibo 
glass have been studied using Ar + and excimer lasers. Non-radiative energy transfer from trivalent 
dysprosium and terbium (donors) to holmium and erbium (acceptors) respectively has been observed 
on the basis of decrease in the life time of the levels and reduction in fluorescence intensity of Dy 3+ 
and Tb 3+ on increasing Ho and Er concentrations. The interaction mechanism of donor and 
acceptor ions is found to be dipole-dipole in both cases. Various parameters such as donor-acceptor 
distances, non-radiative energy transfer efficiency (77) and energy transfer probability (Pda) have 
been computed. Stark splitting have also been marked in several intense transition of the two. 

Keywords. Fluorescence; non-radiative; dipole-dipole; transfer efficiency; Stark components. 
PACSNos 78.66; 78.60 



1. Introduction 

The mechanism of donor-acceptor energy transfer is important in the study of phosphors 
and laser materials. The energy transfer may be through radiative or non-radiative 
processes. The non-radiative energy transfer processes are major factors in determining 
the long wavelength operating limits for the solid state lasers [1]. 

The non-radiative energy transfer has been studied by several authors [2-8] in different 
glasses with different rare-earth ions as dopant. A mechanism for energy transfer from 
donor to acceptor was first suggested by Forster [6] and then extended by Dexter [7]. This 
theory was however found inadequate to explain all concentration regimes of the 
acceptor. Fong and Diestler [8] developed a theory which takes into account the fact that 
for low acceptor concentrations there is usually a one-to-one correspondence between 
donor and acceptor while at higher acceptor concentrations one donor ion may interact 
with two or more neighbouring acceptor ions. 

The non-radiative energy transfer from Dy 3+ to Pr 3 " 1 " in calibo and phosphate glasses 
has been reported by Joshi and coworkers [9, 10]. Similarly energy transfer from Tb 3+ to 
Er 34 " in sodium borate glass has also been studied by Joshi and Pandey [11]. At low 



glass lattice has been studied by Fant et al[15\ and from nr"*" to Er^ in phosphate glass 
by Van Uitert [14] and a dipole-dipole transfer mechanism was reported. 

In this paper we report the fluorescence spectra of Dy, Dy + Ho, Tb and Tb + Er doped 
in calibo glass host. We have been able to resolve the Stark components in the case of the 
intense transitions in Dy 3+ and Tb 3+ . We have also measured the life time of the 
fluorescing levels of the two in the presence and abscence of the acceptor. A significant 
energy transfer from the excited Dy 3+ (Tb 3+ ) to Ho 3+ (Er 3 " 1 ") was noticed from the 
weakening of the fluorescent lines and decrease in life time of Dy 3+ (Tb 3+ ) level in the 
presence of Ho 3+ (Er 3+ ) as acceptor. The mechanism was found to be dominated by 
dipole-dipole interaction. These changes have also been used to estimate the energy 
transfer efficiency (77) and transfer probability (P da ) in the two cases. 



2. Experimental 

The Dy 3+ and Dy 3+ +Ho 3+ , Tb 3+ and Tb 3 + +Er 3+ doped glasses were prepared from the 
spectral grade reagents with 99% stated purity by mixing 20% CaO, 10% Li 2 O 3 and 70% 
B 2 O3. Keeping the Dy 3+ concentration (1 wt%) fixed the concentration of Ho 3+ was 
varied from 0.5 to 3 wt.%. Similarly the concentration of Tb 3+ was held (2 wt.%) fixed 
and Er 3+ concentration was varied from 0.5 to 5 wt.%. The mixture powder was melted in 
platinum crucible heated up to 1000C in an electric furnance and the melt was poured 
into a suitable cast. The moulded glass was grinded and polished carefully to make it 
amenable to spectroscopic studies. 

The density and refractive index were measured and other related physical properties 
were calculated following the procedure outlined in our previous paper [15] and 



Table 1. Various physical parameters of Dy, Dy + Ho, Tb and Tb + Er in calibo glass. 

Donor (Dy 3+ ) and (Tb 3 +) concentration (AT) = 0.42 x 10 20 and 1.02 x 10 20 ions per cc (fixed). 
Donor-donor separation = 28.76 A and 21.04 A for Dy and Tb respectively. 
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Ho 3+ 


Er 3 * 


Cone. (wt%) 


0.0 


0.5 


1.0 


1.5 


2.0 


3.0 


1.0 


2.0 


3.0 


5.0 


Acceptor ion cone. 
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0.45 


0.86 


1.02 


1.56 


2.22 


0.92 


1.56 


2.91 


4.5 
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22.65 


21.40 


18.57 


16.51 


22.15 


18.57 


15.09 


13.22 


separation 






















Acceptor ionic- 
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17.94 


14.45 


13.66 


11.85 


10.54 


14.13 


11.85 


9.63 


8.44 


radius (A) 






















Ref. index 


1.56 


1.53 


1.57 


1.55 


1.52 


1.57 


1.49 


1.50 


1.52 


1.49 


Density 


3.12 


2.95 


3.00 


3.11 


3.10 


2.86 


2.91 


2.87 


3.23 


3.51 


Thickness 


0.15 


0.16 


0.15 


0.14 


0.15 


0.16 


0.25 


0.24 


0.24 


0.26 


Dielectric const. 


2.43 


2.34 


2.46 


2.40 


2.31 


2.46 


2.22 


2.20 


2.21 


2.24 



temperature using Ar + laser lines as exciting line. Ne emission lines from Fe-Ne hollow 
cathode lamp recorded simultaneously were used for wavelength calibration. It is 
interesting to note that the fluorescence intensity increases even up to 6 mol % of the 
Dy 3+ and Tb 3+ in the calibo glass host. 

The life time measurements were made using excimer (XeCl) laser at 308 nm at 10 Hz 
repetition rate and a 5001 model EG&G box car averager with 5 ns gate width. 

3. Results and discussion 

The optical properties of bivalent Dy 3+ and Tb 3+ in glassy and other host matrices have 
been studied extensively [16-18]. From these studies it is well known that Dy 3+ and Tb 3+ 
give broad emissions in the visible region. 

The fluorescence spectra of Dy 3+ , Dy 3+ + Ho 3+ , Tb 3+ and Tb 3+ + Er 3+ in calibo glass 
are shown in figures 1 and 2. Though both Dy 3+ and Tb 3+ ions show fluorescence at 
all the lines emitted by the Ar + laser, the fluorescence yield is maximum at 488.0 nm 
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Figure 2. (A) Fluorescence spectrum of Tb 3+ (2wt%), (B) Tb 3+ (2wt%) + Er 3+ 
(lwt%), (C) Tb 3+ (2wt%) + Er3 + (2wt%) and (D) Tb 3+ (2wt%) + Er 3+ (5wt%) 
doped in calibo glass at 300 K. 
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for Tb 3+ and at 457.9 nm for Dy 3+ . We observed two emission lines for Dy 3+ and 
four lines for Tb 3+ at 25 mW of exciting laser power (see figures 1 and 2). However 
on increasing the laser power to 50 mW some additional lines are seen in the case of 
Tb 3+ (see figure 3), but there is no such effect in Dy 3+ . From the recorded fluore- 
scence profiles at two temperatures (300 and 77 K) in the two cases we notice that at 
liquid nitrogen temperature the peak fluorescence intensity is increased and the bands 
appear sharper. Fluorescence lines and their assignments in the two cases are presented in 
table 2. 

Incorporating Ho 3+ and Er 3+ ions along with Dy 3+ and Tb 3+ respectively (in calibo 
glass) decreases the intensities of the Dy 3+ and the Tb 3+ emission peaks by different 
ratios (see figures 1 and 2). 

3.1 Fluorescence spectra of Dy 3+ and Dv 3+ +Ho 3+ 

The fluorescence spectra of Dy 3+ and Dy 3+ +Ho 3+ in calibo glass excited by the 
457.9 nm line of the Ai + laser are shown in figure 1. The excited Dy 3+ ion from 4 /is/ 2 
state decays rapidly to the 4 F 9 / 2 metastable state. The 4 F 9 /2 level decays by fluorescence 
to the ground 6 Hj multiple! ( 4 F 9 / 2 -^ 6 H^/2 and 4 F 9 / 2 > #15/2) giving emission lines at 
572 nm and 482 nm respectively. The observed line widths are large probably owing 
to inhomogeneous local fields in the glass. One observes the same two lines in 
the fluorescence from other glasses doped with Dy though the observed wavelengths 
in calibo glass are slightly different than for the other glasses (borate, phosphate and 
heavy metal fluoride glasses [19]). It appears that the transition probability from 4 F 9 / 2 
to other components of the ground levels are very small. A look at the fluore- 
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Figure 4. Energy level diagram and energy transfer scheme of Dy 3+ and Ho 3 " 1 " in 
calibo glass. 



as the concentration of Ho increases. No fluorescence lines attributable to Ho are 
observed even up to 800 nm. This latter observation indicates that the Ho level being 
excited through energy transfer from Dy 3 " 1 " decays non-radiatively. The energy separation 
between the 4 F 9 /2 level of Dy 3+ and its ground state (21097cm" 1 ) is almost the same as 
the height of the non-radiative level 5 F 2 (21073 cm" 1 ) of Ho 3+ over its ground state. So 
the energy can be transferred easily from the excited Dy 3+ level to the Ho 3+ level (see 
figure 4). There are a large number of other non-radiative levels in Ho 3+ viz, 3 K%, 5 Ge 
close to the 5 F 2 state while the Ho 3+ states known to decay radiatively have energies 
18650-1 5550 cm" 1 [19-21] above the Ho 3+ ground state. Therefore, the presence of Ho 
ion in the same lattice as the Dy ion depletes the population of 4 F 9 / 2 level of Dy 3+ in a 
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Figure 5. Stark splitting in different energy levels and observed transitions in Dy 3+ 
doped in calibo glass. 



fact that in our study no additional fluorescence lines due to the presence of Ho 
are observed is probably due to 

(i) no direct excitation of the Ho ions as the exciting line is far removed in energy from 

suitable Ho 3+ levels 
(ii) Absorption of the emitted radiation from the Dy 3+ ions is not possible as Ho 3+ 

absorption lines are far removed in energy 
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been observed at considerably high Dy concentrations in various crystalline and glassy 
[20] hosts. The Dy concentration in our case is probably too small to cause any self- 
quenching. 

The two observed fluorescence lines are seen to split into several Stark components 
when the slit width of the spectrometer is reduced and the laser power is increased. The 
spectral peaks observed at 482 nm and at 572 nm show five components each. The 
individual components could be assigned on the basis of the energy level data from Dieke 
[21] and they are shown in figure 5. Their assignments are given in table 3. 



3.2 Fluorescence spectra ofTb 3+ and 

The fluorescence spectra of Tb 3+ doped calibo glass recorded at room and at liquid 
nitrogen temperature show seven emission lines at 491.1, 545.1, 588.1, 624.2, 654.0, 
667.0 and 678. 8 nm corresponding to all the seven components of 1 F J - which is rare. 
They are assigned to 5 Z) 4 > 7 F/ (7 = 6, 5,4, ... ,0) transitions of Tb 3+ and are shown 
in figures 2 and 3. The fluorescence spectrum at liquid nitrogen temperature is sharper 
than the one at room temperature thereby emphasizing the role of the host matrix 
lattice vibrations in causing the line broadening. The fluorescence yield in Tb 3+ is 
maximum for an excitation wavelength of 488.0 nm because the energy of the 5 D 4 
(~ 20462cm" 1 ) level of Tb 3+ above its ground state is closed to the exciting photon 
energy (~ 20491 cm- 1 ). 

Addition of Er 3 " 1 " reduces the intensity of the fluorescence lines of Tb 3+ . At 5 mol % 
of Er 3 " 1 ", the 5 Z) 4 J FI peak of Tb 3+ which appears as a single broad peak at lower 
Er 3+ concetrations split into two well defined peaks. Though 4 F 7 / 2 level of Er 3 " 1 " lie very 
close to the 5 D 4 level of Tb 3+ , it is difficult to say that the extra peak is due to transition 
from 4 F 7 / 2 level of Er 34 " to its lower level. The near equality of the energy difference 
between the 5 D 4 level (~ 20462 cm" 1 ) of Tb 3+ from its ground state to that of the 4 F 7 / 2 
level of the Er 3 " 1 " (see figure 6), permits efficient energy transfer from Tb 3+ to Er 3+ . The 
slight mismatch in the two energy differences may be compensated for by the host lattice 
phonons. The 4 F 7 / 2 level of Er 3 " 1 " is a fluorescing level but it is surprising that it does not 
give any line due to its de-excitation in the spectral region 480-740 nm in calibo host. 
Joshi and Joshi [22] and Joshi et al [23] have made a similar study for energy transfer 
from Eu 3+ to Er 3 " 1 " and from Tb 3+ to Er 3 " 4 " in calibo and sodium phosphate glasses 
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Figure 6. Energy transfer scheme of Tb 3 " 1 " and Er 3 " 1 " in calibo glass. 



observed. It appears that Er 3+ in the 4 F 7 / 2 level relaxes rapidly to the 4 /is/2 level 
(metastable state) which then decays through the transition 4 /is/2 > 4 /is/2 at 1.5 (ira in 
the IR region which is a lasing transition. The energy levels of Tb 3+ and Er 3 " 1 " are shown 
in figure 6. 

The spectral peak observed at 545. 1 nm ascribed to the transition 5 D4> 1 Fs of 
Tb 3+ is very intense and we have succeeded in recording Stark splittings in it by 
increasing the input laser power and by reducing the slit width of the monochromator. 
Five Stark components are seen in this line while for the other 5 D 4 > 7 F 4 transition 
two lines are seen at 583.3 nm and 579.2 nm. Assignments of the observed Stark splitted 
components, have been made on the basis of the known energy levels of the triply 
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Figure 8a. Variation of life time and intensity with acceptor concentration for the 
4 Fg/2 4 /i3/2 transition of Dy 3+ (with Dy 3+ concentration 1% fixed). 
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id absence of acceptor ions Ho 3+ and Er 3+ respectively. The life time of 4 Fg/ 2 level of 
y 3+ were found to be 1.012ms and of 5 D4 level of Tb 3+ 1.485ms in the absence of 
:ceptors. The life time of these levels were found to decrease in the presence of the 
:ceptor. The life time decreases rapidly on increasing the acceptor concentrations. The 
iriation of intensity and life time with acceptor concentration are shown in figures 8a 
id 8b. The absence of any line due to acceptor in the spectra indicates that the energy 
insfer is non radiative. 

The energy transfer probabilities (Pda) and energy transfer efficiencies (r?r) have been 
ilculated employing the formula given by Reisfeld [24] 

Pda = l/T d ((Wfe) - 1) = l/r d ((/do)// d ) ~ 1) = (l/T d - l/T d o) 

rj T = 1 - (W^do) = l ~ (V/do) 

here rjd and ?7do are the donor fluorescence quantum efficiencies in the presence and 
isence of the acceptor respectively (for the same donor concentration) and /d, /do are the 
irresponding fluorescence intensities. Td and r<& are the measured decay time of the 
mor fluorescence in the presence and absence of the acceptor. The measured values are 
ven in table 4. The calculated values for the energy transfer efficiencies and energy 
insfer probabilities are also given in the same table. 

The mechanism responsible for energy transfer from donor to acceptor may be 
entified by plotting energy transfer probability (TV) against the square of the donor + 
:ceptor concentration (Cd + C a ) 2 . This type of plot in the two cases are given in figure 9. 
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of S = 6 indicates that electric dipole-dipole coupling is the predominant interaction 
mechanism for energy transfer. At very large concentration of acceptor a deviation from 
the linearity is marked. Such a deviation may be due to some other mechanism which has 
started taking place in the glass with the larger acceptor concentration. This is also 
supported by an abrupt decrease in fluorescence intensity of Tb 3+ when the concentration 
of Er 34 - is 5%. 



5. Conclusion 

Fluorescence spectra of Dy 3+ , Dy 3+ + Ho 3+ , Tb 3+ and Tb 3+ + Er 3+ have been studied at 
different concentration of Ho 3+ and Er 34 ". From the variation in intensity of the spectral 
lines it is concluded that Dy 3+ and Tb 3+ are good sensitizers for Ho 3+ and Er 34 " ions in 
calibo glass hosts. The lifetime of the 4 F 9 / 2 level of Dy 3+ and 5 >4 level of Tb 3+ in 
presence and absence of acceptor Ho 3+ and Er 3+ respectively have been measured. From 
this it is concluded that a non-radiative energy transfer is taking place and the mechanism 
is electric dipole-dipole. Stark splitting have been marked in some of the strong 
transitions of the two. 
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Abstract. Two-photon optogalvanic transitions in Ar glow discharge with Nd : YAG laser pumped 
dye laser excitation in the frequency range 13520-1 6520 cm" 1 has been studied using linear 
and circular polarization. The intensities of two-photon optogalvanic transitions are very sensitive 
to changes in the incident laser power which is not the case with one-photon transitions. Intensity 
ratio for circular and linear polarized light for two photon transitions 6j'[l/2]l < 4y[3/2]2, 
6s'[l/2]0<-4s[3/2] 2, and 5d[l/2]0 <- 4j[3/2] 2,5d[l/2] l <- 4s[3/2]2 are quite different 
from the other two-photon transitions. This has been explained as due to near one photon resonance 
of 4 J p'[3/2]l level for the first pair and 4p'[l/2]l for the second pair of transitions. The ratio of 
optogalvanic intensity for circular to linear polarized light has been theoretically estimated and 
compared with the observed results. 

Keywords. Optogalvanic spectrum; two-photon transitions; laser spectroscopy. 
PACSNos 52.80; 42.60 

1. Introduction 

Optogalvanic spectroscopy of glow discharges is based on the detection of small changes 
in the plasma-impedance, with high signal to noise ratio. This change of impedance is 
caused by resonant absorption of monochromatic radiation from the tunable laser, by 
atomic or molecular species in the plasma. Hollow cathode discharges have been found 
particularly suitable for "such studies on the excited electronic states and rare gas filled 
commercial hollow cathode lamps have also been used for the wavelength calibration of 
dye lasers [1-3]. Extensive studies of neon optogalvanic spectra have led to a great deal 
of information on its one-photon and two-photon transitions [2, 3]. In contrast, there are 
no reports on the two-photon optogalvanic transitions in argon in spite of several studies 
of its one-photon optogalvanic spectra [4-7]. There are also no systematic investigations 
of the extreme red end of the Ar optogalvanic spectrum apart from identification of about 
ten one-photon lines spread over the 695-795 nm region [4]. We report here, for the first 
time, two-photon optogalvanic transitions in Ar with Nd : YAG laser pumped dye laser 
excitation in the wavelength transitions region 605-740 nm (16520-13520 cm" 1 ). The 



(model YG681C) pumped a Lambda Physik dye laser (model LPD3002) at 10 Hz with 
energy per pulse 250 mJ. Three different laser dyes in methanol were used: LDS-730 (740- 
720 nm, region I), LDS-698 (720-670 nm, region II) and DCM (670-605 nm, region IE). 
The spectral bandwidth of the dye laser was 0.2 cm" 1 with energies (at the peak of the dye 
laser output) per pulse (5 ns duration) of 35 mJ in region I, 70 mJ in region II and 45 mJ in 
region IH Approximately 20% of the dye laser radiation was coaxially directed into a 
Hamamatsu Ar-hollow cathode lamp (model L233-Fe) which was connected to a regular 
d.c. power supply (Keithley model 246) through a current limiting resistance of 2.5 kQ. 
The voltage across the lamp could be set between 180 and 200 volts for normal glow 
discharge and the optogalvanic voltage change of the discharge was monitored across a 
2.5 k O resistance using a coupling capacitor of 5 u,F. The laser beam spot was smaller than 
the cathode diameter of the lamp and the unfocused laser beam was carefully aligned so as 
not to fall on the metal surface. Our attempts to focus the laser beam (even weakly), in the 
glow discharge, led to pronounced instabilities in the optogalvanic signal. 

The optogalvanic signal were monitored on an oscilloscope (Tektronix model 2465A) 
and simultaneously processed by a boxcar averager (model SR250); both units were 
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Figure 1. Oscilloscope traces of optogalvanic signals in argon corresponding to 
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Figure 2. Optogalvanic spectra of Ar with a boxcar gate delay of (a) 13 us (linear 
polarization) (b) 2 us (linear polarization) and (c) 2 us (circular polarization). 



triggered by a photodiode receiving a small fraction of the incident dye laser pulse. A 
general survey of the temporal profile of the signals at several wavelengths showed that 
they could be grouped into two types as shown in oscilloscope traces in figure 1. The first 
group of transitions exhibit a negative voltage profile with a minimum at 1 us from the 
excitation laser pulse and asymptotically increase to after about 20 u,s (figure Ic, d). The 
second group of optogalvanic transitions exhibit a composite voltage profile which starts 
with the negative voltage signal, becomes minimum after 2^is then increases to a 
maximum positive voltage after 10 to 15 \JLS and asymptotically decreases to at the end 
of 50 [is (figure la, b). The boxcar processed optogalvanic signals were recorded on a 
chart recorder with the dye laser grating scan rate of 0.3 nm/min. Figure 2 shows the 
optogalvanic spectra in the frequency region 13950-14400 cm"" 1 (7 17-694 nm) recorded 
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Figure 3. Optogalvanic spectrum of argon with a boxcar gate delay of 4 \is (linear 
polarization). The numbers within brackets correspond to serial number of transitions 
as in table 1 and (*) denote two-photon lines. 

transitions are greatly diminished for the longer boxcar delay (figure 2a) and the 
Optogalvanic signals for the second group become positive in polarity. It was found that 
the maximum number of Ar lines (all with the negative voltage signals) in any dye laser 
region were simultaneously recorded by setting the boxcar gate delay between 2 u,s and 
4u.s. The spectrum shown in figure 3 in the frequency region 14900-1 5400 cm" 1 (650- 
670 nm) is a part of the Optogalvanic spectra in the DCM dye region (region III) recorded 
with the gate delay of 4 u.s. 



3. Spectroscopic assignments 

The frequencies of the observed spectral lines were determined by fitting a second degree 
polynomial to the measured line positions with a maximum uncertainty of 1 cm" 1 , these 
are given in the first column of table 1. The relative Optogalvanic intensities (estimated 
from peak height for linearly polarized light) are given in the second column but these are 
not normalized for the spectral power variations of the dye laser output and in the case of 
overlapping transitions the relative intensities in the two dye regions are also shown 
(within bracket). The upper and lower energy levels of the transition, in Racah notation 
[8], are given in the third and fourth columns respectively and the transition type (one or 
two-photon) is shown in the fifth column. In region n, the relative intensities of two- 
photon transitions were also determined for circularly polarized light and the ratio (0) of 
circular to linear polarization intensities is included in the fifth column along with the 
corresponding uncertainty. The spectroscopic assignments are quite straight forward in 
terms of the known energy levels of argon [9] and all, except ten of 142, observed lines 



Assignment 



Frequency 
(cm- 1 ) 


rs.ciau.ve 

intensity 


Upper state 


Lower state 


llCtllMUUll 

type (fi) 


Region I (Dye 


LDS-730) 








13523.2 


03 


643/2] 1 


4p[3/2]2 


OP 


13539.9 


18 


4p'[3/2]2 


443/2] 1 


OP 






743/2] 1 


4s' [1/2]! 


TP 


13561.4 


10 


4rf[7/2]4 


4p[5/2]3 


OP 


13595.3 


12 


4d[l/2}3 


4p[5/2}2 


OP 


13610.5 


04 


4d[3/2]l 


4p[3/2}2 


OP 


13664.8 


07 


6s'[l/2]l 


4p'[l/2]l 


OP 


13673.2 


04 


643/2] 1 


4p[3/2]l 


OP 


13705.4 


14 


6s'[l/2]l 


443/2] 1 


TP 


13722.1 


05 


4d'[3/2}\ 


V[3/2]2 


OP 


13745.8 


35 


4p'[l/2]l 


4^[3/2]l 


OP 


13749.9 


17 


4d[7/2}3 


4p[5/2]3 


OP 


13757.1 


03 


743/2] 1 


V[l/2]0 


OP 


13760.5 


05 


4d[3/2]l 


4p[3/2]l 


OP 


13804.9 


06 


4d'(5/2}3 


4s[3/2}2 


TP 


13827.5 


09 


4d[5/2}2 


4p[5/2}2 


OP 


13872.0 


17(b) 


6/[l/2]l 


V[3/2]2 


OP 


13880.0 


05 


4d'[3/2}\ 


V[3/2]l 


OP 


13930.0 


05(b) 





unassigned 




13944.4 


05 


743/2]2 


4s / [l/2]0 


TP 


13948.3 


07 


4d[5/2}3 


4p(5/2}2 


OP 


Region H (Dye-LDS-698) 


13976.6 


70 


6J[l/2]Q 


4s[3/2}2 


TP (0.9 0.1) 


13988.3 


60 


V[3/2]l 


4^[3/2]2 


OP 


14008.8 


80 


6j'[l/2] l 


4j[3/2]2 


TP (1.3 0.1) 


14079.2 


06 


5d'[5/2}3 


4/[l/2]l 


TP (2.0 0.1) 


14091.6 


11 


5rf[l/2]l 


443/2] 1 


TP (0.7 0.2) 


14104.3 


04 


4d[5/2}3 


4p[5/2}3 


OP 


14128.2 


06 




unassigned 




14145.9 


100 


4p'(3/2}2 


4^[3/2]2 


OP 


14168.4 


07 


5d[3/2}2 


443/2]l 


TP (2.1 0.2) 


14204.5 


06 


5d[J/2]3 


443/2]l 


TP (1.9 0.2) 


14222.3 


03 


4rf[5/2]3 


4p[5/2}3 


OP 


14231.8 


02 


4d[3/2]l 


4p[5/2}2 


OP 


14266.5 


04 


5d'[5/2]2 


443/2]l 


TP (2.1 0.1) 


14290.5 


02 


5d'[5/2}3 


443/2]l 


TP (2.0 0.1) 
HP 



intensity 



Upper state 



Lower state 



type (Q) 



14324.9 


100 


5d[l/2}0 


4 S (3/2}2 


TP (1.1 0.1) 


14352.7 


65 


4p'[l/2]l 


4^[3/2] 2 


OP 


14365.0 


15 


7j[3/2]l 


4s[3/2]l 


TP (0.9 0.2) 


14382.0 


06 


5d[3/2]l 


4s[3/2]l 


TP (0.8 0.2) 


14393.7 


65 


5d[l/2]l 


4^[3/2]2 


TP (1.5 0.1) 


14410.0 


07 


5d'(3/2}2 


4s'(l/2}Q 


TP (2.0 0.2) 


14445.9 


13 


5d[l/2}4 


4^[3/2]2 


TP (2.1 0.2) 


14471.7 


20 


5d[3/2}2 


4s[3/2] 2 


TP (0.7 0.2) 


14475.7 


05 


5d'[5/2}2 


4^'[1/2]0 


TP (1.9 0.2) 


14508.7 


03 


5d{l/2}3 


4^[3/2]2 


TP (2.0 0.2) 


14515.3 






(unassigned) 




14550.5 


03 


4d[l/2]l 


4p[l/2]l 


OP 


14592.8 


15 


5d[3/2}2 


4p'[l/2]l 


OP 


14592.8 


04 


5d[5/2}3 


45[3/2]2 


TP (2.0 0.2) 


14648.7 


09 


ls[3/2}2 


4j[3/2]2 


TP (0.7 0.2) 


14785.7 


10 


5d[5/2}2 


4p'[l/2]l 


OP 


14803.1 






(unassigned) 




Region ffl 


(Dye-DCM) 








14924.5 


25 


4d'[3/2}l 


4p[3/2]l 


OP 






6s'[l/2]l 


4p[3/2]2 


OP 


14944.5 


08 


ls[3/2}2 


4p'[l/2]l 


OP 


14956.2 


20 


5d[3/2}2 


4p'[3/2]l 


OP 


14972.5 


40 


V[l/2]0 


4*[3/2]l 


OP 


14983.5 


18 


4d'[3/2}2 


4p[5/2}2 


OP 






ls[3/2]l 


4p'[l/2]l 


OP 


15002.0 


35 


4d'[5/2]2 


4p(5/2}2 


OP** 


15010.1 


04 


6a / [l/2]0 


4P[3/2]1 


OP 


15018.1 


08 


5d[3/2]l 


4p'[l/2]l 


OP 


15039.7 


10 (b) 


5d[5/2]3 


4p'[3/2]2 


OP 


15065.5 


09 


7s / [l/2]l 


4s[3/2]l 


TP 


15074.7 


09 


6s'[l/2]l 


4p[3/2]l 


OP 


15093.2 


08 


8j[3/2]l 


4j[3/2]l 


TP 


15114.5 


30 


5d'[3/2]2 


4^[3/2]2 


TP 


15136.2 


70 


4d'[5/2]3 


4p[5/2}2 


OP 


15137.8 


75 


4d'[3/2]2 


4p(5/2}3 


OP 


15150.2 


28 


4d[5/2}2 


V[3/2]l 


OP 


15155.7 


60 


4d'[5/2]2 


4p[5/2]3 


OP** 


15161.2 


06 


6rf[l/2]l 


4s[3/2]2 


TP 


15181.8 


06 


6d[l/2]0 


4s[3/2]2 


TP 


15190.1 


06 


ls[3/2}l 


4p / f3/2l2 


OP 



p 


Rel stive 






Transit'on 


(cm- 1 ) 


intensity 


Upper state 


Lower state 


type (ft) 


15213.7 


09 


8j[3/2]l 


V[l/2]0 


OP 


15224.8 


06 


5d[3/2]l 


V[3/2]2 


OP 


15254.0 


20 


6d[7/2}4 


4^[3/2]2 


TP 


15273.6 


07 


10J[5/2]2 


4s > [\/2\\ 


TP 


15279.2 


07 


6d'[5/2}2 


4s / [l/2]0 


TP 


15290.4 


90 


4d'[5/2]3 


4p[5/2}3 


op** 


15308.1 


07 


7s[3/2}2 


4p'[3/2}l 


OP 


15342.2 


25 


4d[5/2]2 


4j?'[l/2]l 


OP 


15347.0 


12 


7s[3/2}2 


V[3/2]l 


OP 


15379.1 


14 


8^[3/2]2 


4^[3/2]2 


TP 


15382.5 


12 


5d[3/2]l 


V[3/2]l 


OP 


15394.4 


28 


4d'[3/2]l 


4^[5/2]2 


OP 


15426.1 


05 


7s[3/2]l 


4p[l/2]0 


OP 


15449.2 


03 


7d[7/2}3 


4s[3/2]l 


TP 


15460.8 


08 


5d[3/2}l 


4p[l/2]0 


OP 


15491.7 


03 (b) 




unassigned 




15515.9 


03 


9s[3/2}l 


443/2] 1 


TP 


15544.3 


25 


6j?'[l/2]l 


4p[5/2}2 


OP 


15580.9 


100 


6s[3/2}2 


4p[l/2]l 


OP 


15624.8 


02 (b) 




unassigned 




15641.7 


03 (b) 




unassigned 




15657.3 


75 


6s[3/2]l 


4p[l/2]l 


OP 


15695.1 


50 


5d[l/2]"l 


4p[3/2]2 


OP 


15706.4 


. 30 


5d[l/2]0 


4p[3/2]l 


OP 


15729.8 


08 


7d[3/2]2 


4j[3/2]2 


TP 


15732.8 


14 


7d[7/2}4 


4j[3/2]2 


TP 


15744.8 


35 


4d[3/2]l 


4p[l/2]l 


OP 


15785.1 


03 


7d[5/2]3 


4s[3/2}2 


TP 


15790.2 


03 


10j[3/2]l 


4s[3/2]l 


TP 


15813.6 


13 


9s[3/2]2 


4s[3/2]2 


TP 


15832.6 


07 


7d[l/2]l 


V[l/2]0 


OP 


15845.6 


50 


5d[l/2]l 


4p[3/2]l 


OP 


15849.8 


70 


5d[3/2}2 


4p[3/2]2 


OP 


15877.0 


25 


5d'[3/2]2 


4p'[l/2]l 


OP 


15922.6 


40 


5d[7/2}3 


4p[3/2]2 


OP 


15961.7 


04 


6d'[7/2}2 


4 5 [3/2]2 


. Tp 


15972.1 


15 


6d[l/2]l 


4p'[l/2]l 


OP 


16000.7 


25 


5d[3/2}2 


4p(3/2]l 


OP 


16009.9 


12 


8d[l/2]0 


4s[3/2]2 


TP 


16013.0 


12 


6d[l/2l0 


4p'[l/2]l 


OP 



or 1 ) 


intensity 


Upper state 


Lower state 


type (J2) 


>045.4 


35 


5d[5/2}2 


4p[3/2]2 


OP 


i083.5 


60 


5rf'[3/2]2 


4p'[3/2]2 


OP 


J092.9 


75 


5d[5/2}3 


4p[3/2]2 


OP 






10s[3/2]2 


4^[3/2]2 


TP 


5127.2 


03 




unassigned 




5165.5 


25 (b) 




unassigned 




5178.8 


13 


6d[l/2]l 


4p'[3/2]2 


OP 


5195.5 


40 


5d[5/2}2 


4p[3/2]l 


OP 


5202.7 


55 


ls[3/2}2 


4p[3/2}2 


OP 






lld[7/2]3 


4j[3/2]l 


TP 


5216.4 


40 


5d'[5/2}2 


4p'(3/2}2 


OP 


5242.4 


50 


7s[3/2]l 


4p(3/2}2 


OP 


5267.8 


30 


5J'[5/2]3 


4p'{3/2}2 


OP 


5276.7 


17 


5d[3/2]l 


4p(3/2]2 


OP 


5312.3 


14 


6d[3/2]2 


4p'[l/2]l 


OP 


5315.5 


28 


5rf[l/2]l 


4p[5/2]2 


OP 


6318.7 


14 


5cf'[3/2]l 


4p'[l/2]l 


OP 


6329.9 


13 


6d[5/2}2 


4p'[l/2]l 


OP 


6336.3 


14 


6d[l/2]l 


4p'[3/2]l 


OP 


6353.1 


13 


7s[3/2]l 


4p[3/2]l 


OP 


6376.3 


50 


6d[l/2]0 


V[3/2]l 


OP 


6386.0 


45 


7j / [l/2]l 


4p'[l/2]l 


OP 


6391.2 


40 


7j[3/2]l 


4p[3/2]l 


OP 


6406.2 


05 


85[3/2]2 


4p'[l/2]l 


OP 


6413.3 


30 


6d[l/2]l 


4p(l/2}G 


OP 


6426.2 


10 


5d[3/2]l 


4p[3/2]l 


OP 


6439.1 


10 


8s[3/2]l 


4p'[l/2]l 


OP** 


6453.0 


15 (b) 




unassigned 




6483.9 


12 


6d[7/2]3 


V[3/2]2 


OP 


6498.5 


25 


4d'(3/2}2 


4p[l/2]l 


OP** 


6516.0 


18 


4d'[5/2}2 


4p[l/2]l 


OP 



b): Broad feature; **: Revised assignment from those in Ref [6]. 



The metastable 4s[3/2]2 level accounts for more than half of the 44 two-photon 
esonant transitions whereas 14 of them originate in 4s[3/2]l, four in 4s / [l/2]0 and 
inly three in 4s t [l/2} 1 level. A partial energy level diagram of argon is given in figure 4 
o explain the origin of some very intense spectral lines. The intensities of two- 
ihoton optogalvanic transitions are very sensitive to changes in incident laser power, 
inlike the one ohoton transition. The srronpp.st twn-nhntnn lin^c nr^ rthooruH ot 



Table 2. Intensity ratio (fi) circular to linear polarized radiation of the two-photon 
transitions. 



ition 
jncy 

') 


Assignment 


Observed 
(ft) 


Calculated 
(0) 


Upper state 


Lower state 


i.6 


6s'[l/2]0 


4s[3/2]2 


(0.9 0.1) 


1.5 


1.8 


6s'[l/2]l 


4s[3/2]2 


(1.3 0.1) 


1.5 


.2 


5d'{5/2}3 


4s f [l/2]l 


(2.0 0.1) 


1.5 


..6 


5d[l/2]l 


4^3/2] 1 


(0.7 0.2) 


1.0 


!.4 


5d[3/2}2 


4^3/2]l 


(2.1 0.2) 


1.5 


L5 


5d[l/2}3 


45[3/2] 1 


(1.9 0.2) 


1.5 


i.5 


5d[5/2}2 


4j[3/2]l 


(2.1 0.1) 


1.5 


).5 


5d'[5/2}3 


443/2]l 


(2.0 0.1) 


1.5 


L9 


5d[l/2]0 


4s[3/2]2 


(1.1 0.1) 


1.5 


i.O 


743/2]l 


443/2]l 


(0.9 0.2) 


1.0 


1.0 


5d[3/2}l 


443/2] 1 


(0.8 0.2) 


1.0 


i.7 


541/2] 1 


4s[3/2]2 


(1.50.1) 


1.5 


).0 


5d'[3/2}2 


4y[l/2]0 


(2.0 0.2) 


1.5 


i.9 


5d[l/2}4 


45[3/2]2 


(2.1 0.2) 


1.5 


..7 


5d[3/2}2 


4y[3/2]2 


(0.7 0.2) 


0.6 


j.7 


5d'[5/2}2 


4s / [l/2]0 


(1.9 0.2) 


1.5 


5.7 


5d[7/2}3 


4s[3/2]2 


(2.0 0.2) 


1.5 


>.8 


5d[5/2}3 


45[3/2]2 


(2.0 0.2) 


1.5 


1.7 


7s[3/2}2 


4r[3/2]2 


(0.7 0.2) 


0.6 



ks[3/2]2 level. The relative intensities of these transitions under linear and cir- 
polarizations are shown in figures 2b and 2c respectively and discussed in the next 
:>n. 

ily ten of the one-photon transitions have been observed in the region II (table 1) of 
h three have metastable states as their lower energy levels. They originate in 
'2]2 level and are very strong, probably power broadened due to high laser intensity 
figure 2). The transitions originating in 4p[5/2]2, 4p[5/2]3 and 4p[l/2]l give rise to 
: optogalvanic signals. All the one-photon optogalvanic transitions in regions I and II 
nany of them in region in have been observed for the first time; for some of the 
itions in region HI reported earlier [6] we have noticed printing or assignment errors 
ked ** in table 1). 



So/Miccisin anrl 



range of 10-20 ns. Temporal profile of the optogalvanic voltage signal resulting from a 
transition between two excited states is given by [11]: 

AV = VQ(HI n2)[fl2exp( t/T-i} a\ exp( t/T\}\, (1) 

where ai, a\(ai > a\] are characteristic phenomenological constants, TI, T\ the relaxa- 
tion times and HI, n\ the atomic populations of the upper and lower energy levels 
respectively. VQ is a constant for a particular transition and depends on its optical 
absorption cross section as well as the discharge and circuit parameters. 

In temporal profile of the optogalvanic signals immediately after the laser excitation 
pulse (t = 0) the voltage is negative because of the dominance of the first term within the 
square brackets of eq. (1). This is exemplified by all the transitions shown in figure 2b. 
For the optogalvanic transitions with 4s[3/2]2 as the lower state and 4p'[3/2]2, 4/?'[3/2]l, 
4p'[l/2]l as upper states, the signals are positive for long time delay (figure 2a) but 
negative for short delay (figure 2b) after the exciting laser pulse. The composite negative 
and positive optogalvanic signals appear because of T% <T\, so that the first term 
in eq. (1) dominates at short delays and the second term becomes more dominant for 
longer delays after the exciting laser pulse. This implies that the upper states of these 
transitions have relaxation times comparable to those of the lower states in accordance 
with eq. (1). 

The observed ratio fi of two-photon intensities with circular and linear polarization 
can be broadly grouped into two categories. For transitions with equal J for the upper 
and lower states fi w 0.7 and for those with unequal / values, it is close to 2. Four of 
the strongest two-photon lines, however, do not show this polarization behaviour (see 
table 1). 

Grynberg et al [12] have carried out a detailed analysis of optogalvanic intensity, 
arising from the scalar and the quadrupolar components of the two-photon transition 
tensor for Doppler-free two-photon spectrum of neon. They conclude that the experi- 
mental observations on polarization characteristics do not agree with theoretical predic- 
tion unless one makes use of accurate wavefunctions for the initial, final and various 
intermediate states. We also qualitatively account for the observed ratio (Q,) in the 
-following way. We have calculated the optogalvanic intensity ratios of circular to linear 
polarized light (1 = circular /linear) by evaluating the contributions of scalar and 
quadrupolar components using the two-photon matrix given by Grynberg et al [12]; 
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and hence 
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where 



1 ^(circular) / 1 o(Iinear) 

1/5C 
-2/15C" 

B = \(eJ e \\Q\\sJ 8 }\ 2 , 
C=\(eJ e \\Q 2 \\8J 8 }\ 2 - 

The following derivations of 0, are valid only when there is only one intermediate real 
level with energy close to that of one photon. 

Case 1. When J e is different from J 8 , then |{e/||>||g./}| 2 = so that 

1 /C 

1.5. (8) 



2. When J e is equal to J g and greater than or equal to 1, then 

R= \(eJ e \\Q\\gJ S }\ 2 = B 
15^!2 22 1/5C' 



(9) 



where R is the ratio of scalar to quadrupolar components, the value of R has been 
calculated from the Racah coefficients according to the following relation: 



R(g,eJ r } = 
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(ii) 



The observed values of fi for transitions at 14471.7cm l and 14648.7cm l as shown in 
table 2 are in good agreement with the calculated value of 0.6 corresponding to transitions 
between two levels with J' = J" = 2. Similarly the observed values of fi for transitions at 
14365.0cm" 1 and 14382.0cm" 1 agree well with the calculated value of 1.0 corresponding 
to transitions with J 1 = J" = 1. However, the calculated value of 1 for transitions corres- 
ponding to J' ^ J" do not agree with the observations. except for that at 14393.7 cm" 1 with 
il = 1.5 0.1 compared to the theoretical value of 1.5. Thus it is found that a single 
intermediate level approximation is not very successful in the calculation of O for two 
photon transition intensities between levels with unequal J values. 

The enormously large intensity of two pairs of two-photon transitions shown in 
figures 2b and 2c must be due to the near resonance intermediate state (figure 4). The 
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Figure 4. Partial energy level diagram of Ar showing (i) the 3p 5 4s metastable states, 
(ii) the strongest two-photon (heavy arrows) and one-photon (light arrows) transitions 



tions around 14353cm l (figure 4). Biraben et al [13] have also made parallel observa- 
tions in optogalvanic spectrum of neon. The effects of linear and circular polarization on 
these two-photon optogalvanic signals are shown in figures 2b and 2c respectively. The 
one-photon lines located between the two pairs of two photon transition may be used as 
calibration marks for the change in relative intensities of the two-photon lines. The two- 
photon intensity ratio fi (circular/linear) is found to be 1.5, 1.1, 1.3 and 0.9 for transitions 
to 5rf[l/2]l, 5d[l/2]0, 6s'[l/2]l and 6,v'[l/2]0 upper levels respectively from the 
metastable 4.s[3/2] 2 level. This observation does not follow the polarization behaviour 
found for other two-photon transitions with unequal J for upper and lower states (see 
preceding paragraph). The near resonance of the intermediate states and the important 
role of collision induced ionization [14] in the generation of optogalvanic signals may be 
responsible for unusual polarization characteristics of these transitions. 

5. Conclusion 

One-photon transitions exhibit a composite temporal profile with both negative and 
positive voltage signal whereas the two-photon transition profiles have negative voltage 
signal only. Linear and circular polarization studies have been carried out to distinguish 
different types of two-photon transitions. 

A near resonant observation has been made for the first time where intermediate level 
in the two-photon transition happens to be a real level. This gives rise to intensity changes 
in the circular and linear polarizations as compared to two-photon transitions with 
nonresonant intermediate levels. 
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Abstract. An electrostatic charge state selector has been constructed for charge transfer studies in 
ion atom collisions. Its design and performance have been discussed illustrating with examples of 
some data taken using heavy ion beams from the pelletron accelerator. Expressions for the deter- 
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operation and also the line shapes have been discussed analytically in detail. 
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1. Introduction 

In ion-atom collisions several single or multi-electron transfer processes occur between the 
projectile and the target atoms resulting invariably in a change of the charge carried by the 
incident ion. It is often required to associate the observable event in the collision process to 
the corresponding charge state of the ion after the collision has taken place. Thus analysis 
of the charge state of the ion after the collision becomes necessary. The final charge state 
information on the ion also provides vital information in some special cases. In the case of 
energetic heavy ions channeled through a thin crystal, for example, the 'frozen' fraction of 
the incident charge state is an important parameter in deriving capture cross sections (see 
e.g. ref. [1]). In recent years our group has undertaken a series of investigations concerning 
various aspects of ion-atom collisions in a solid target. A charge state selector was required 
to resolve the charge states of the ion after the collision process. We describe here the 
construction of a simple electrostatic charge state selector and illustrate its performance 
with examples of some of the data taken using BARC-TTJFR pelletron accelerator. 
Analytical derivations of the possible line-shapes of the charge spectrum and procedure for 



voltage up to 50 kV could be applied. One of the plates was arranged at an inclination 
with respect to the other in order to avoid restrictions imposed by the parallel plate 
geometry thus effectively achieving higher charge dispersion than in a parallel plates 
arrangement. The plates were mounted in such a way that a positively charged particle is 
deflected in the vertical direction. The deflector electrodes were made from highly 
polished 74.5 cm long, 30 mm wide, and 0.5 cm thick SS plates. The distance between the 
plates was maintained to be 16mm and 42mm at the entrance and exit, respectively, by 
means of three sets of appropriately machined ceramic insulators. A drift space of 5 1 cm, 
which also accommodated the pumping system, was kept between the exit of the 
deflector plates and the Faraday cup. The vertically deflected ions were collected in a 
electron suppressed Faraday cup placed just behind a 2 mm high and 10 mm wide slit kept 
at a vertical distance of 45 mm from the undetected beam position. A slit 1 .5 mm wide 
and 10mm high was placed at the entrance of the charge state selector to restrict the 
acceptance angle of the beam. A line drawing of the charge selector assembly is shown in 
figure 1. The high voltage on the plates was provided by a shielded cable which was 
connected via specially made high vacuum teflon feed-throughs (figure 2). The Faraday 
cup, behind the slits, was directly fitted to a standard vacuum feed-through which was 
also electrically isolated. Two externally programmable high voltage supplies were 
employed for providing the required voltages to the plates. 
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3. Control and data collection 

3.1 Modes of operation of the charge state selector 

The charge spectra can be obtained either by recording the spectra at a fixed plate \ultaj?r 
in conjunction with a position sensitive particle detector at the exit of the charge, selivtui 
(to be incorporated in near future) or, as is clone in the present case, by recording ihr 
charge transmitted through a narrow vertical slit kept at a fixed distance, as a luwliun ui 
the plate voltage. Let us call the first mode of operation, PSD mode and the present urn*, 
voltage scanning (VS) mode. In VS mode, the charge is collected via a speeiaih designed 
Faraday cup, kept just behind the slit. The charge collected by the Faraday vup is 
measured using an electrometer whose output was digiti/ed. For incident charge 
normalization, the electrons emitted by the rotating wire of a beam profile moniuu 
(BPM) mounted before the scattering chamber were collected and measured in similar 
way by another electrometer. This arrangement was particularly suitable while usinj.* vei> 
small currents which are inevitable in the case of fully stripped ions. 



3.2 Hardware 

In order to vary the voltage on the plates of the charge selector in finite steps and otllevt 
the data at each step, a PC based system was designed. Figure 3 slums the How chart and 
the functional block diagram of the control and data collection .system. 

The interface (I/F) board consists of address, data and control signal huttei It's. 
(2 x 74 ALS 244, 74 ALS 245), an address decoder (74 AI .S 1 3K ), magnitude coiujMiator 
74LS6881, two 8255 programmable peripheral I/O devices, two KIM count duvui 
counter ICs, a 10 bit D/A converter (AD56U) along with a last settling time t *HM 
operational amplifier (LF396) along with a sample and hold amplifier tl.I- WK. Ttu* 
buffered address lines, BA5 to BA9, are used to select the group address < MX* Hex in uw 
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Figure 3. Row chart and the functional block diagram of the control and data 
collection system of the charge state selector. 



address lines BAG and BA1. Each of the two 8255 ICs on the I/F board has three ports 
(A, B and C). One of the 8255 (CS1) is used to operate the D/A converter (8 bits of 
port A and 2 bits of port B together form 10 bits to drive the DAC). The output from the 
DAC is given to the operational amplifier LF 356, which gives a -f-ve output in the range 
of (0-5 V) or (0-10 V) selectable by a dip switch. An inverter is provided for ve analog 
voltages of the same range. A sample and hold (S/H) amplifier serves to generate a 
steady DC signal for driving the programmable HV supplies. The functions of the S/H 
amplifier are: (i) to isolate the output signal from the input signal while the input bits of 



the high voltage supply from software mistakes and hardware faults which would lead to 
full 50 kV on each plate causing sparks. We have provided one more 8255 (CSO) for 
control and operation of other devices (like stepper motors etc.) in future. 

In order to collect data from various devices like Faraday cups, detectors, current 
integrators, timers etc., six 16 bit count down counters are provided using two 8254 ICs 
(CS2 and CSS). As soon as the plate voltage is set, the counters are loaded with maximum 
count length (max value of 16 bits) and are enabled. These counts are down counted. The 
data is fed at the i/p of each counter in the form of serial TTL pulses. After a preset time 
or counts is reached, the counters are disabled, the data from the counters are read to 
obtain the actual counts and stored in files. All the data IN and OUT are transferred using 
the buffered data lines BDO-BD7. The read and write instructions from the PC are 
given by the BIOR and BIOW lines through the decoder. Simultaneously, the plate 
voltage is incremented to the next value and the cycle is repeated. The counters are 
enabled or disabled using the PCI and PCO lines of the port C of the first 8255 (CS1) 
device. After the data collection at the highest preset value of the plate voltage, the 
complete data is transferred to a previously named file in the PC. Also the data from 
counters 1 and 3 are displayed on the monitor. It is also possible to interactively 
normalize the data of counter 3 to the corresponding data counter 1 and display the 
normalized data. The whole set of data collection with same parameter can be repeated 
interactively by a single key. 

3.3 Software 

A user friendly menu driven interactive software using 'Basic' was written. This software 
provides several options, which are briefly described below: 

Scan: Using this option the plate voltages (both positive and negative) can be set to vary 
simultaneously with a common control between two limits in finite steps. Six set of 
charge data collected on various slits, FC etc. or particles detected in the counters can be 
recorded at each step for a given time interval or for a preset number of incident particles. 
Once the selected voltage scan is completed, the number of contents of the two counters 
are displayed as a function of the applied voltage. Simultaneously, the data from all 
counters are stored in a file. The data of one of the counters could be normalized with 
respect to the other counter at each voltage step and can be viewed. The scan with same 
input parameters can be repeated several times either automatically or interactively. 

Display: This option allows the data of two files stored earlier to be displayed. 

View: In this option the data from all counters recorded earlier in a file along with various 
parameters used are displayed on the monitor screen of the PC. 

Files: Various files existing in a directory can be viewed using this option. 



beam of charge q and energy E in the plates of length / and a drift region /drift can be 
expressed as 

d e / dri ft 
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2 d e do [ I * - o ' J V*o/ J 

(1) 

where V is the applied voltage difference between the plates and d$ and d e are the 
distances between the plates of the charge selector at the entrance and exit, respectively. 
yi nc is the transverse displacement of the charged particle from the optic axis and a\ nc is 
the angle of incidence with respect to the optic axis. Let us introduce the constants C\ and 
Ci so that s can be expressed as 

s = }>inc + Ci tan a inc + C 2 ( ) . (2) 



C\ and Ci are functions of the geometrical parameters of the charge state selector only. 
The deflection s can be obtained directly in mm when one expresses /, d e> do, and /drift in 
mm, the projectile energy E in MeV, the voltage V in kilo- Volts and the projectile charge 
q in eletronic unit. 

As already mentioned, the charge state of the incident ion may change during the 
collision with the target atoms and a distribution of charge states is produced after the 
collision. In this section we derive the fractions of post-collision charge states which enter 
the selector and also their spectral line shapes. 

Typical charge spectra normalized to the number of ions incident on the target are shown 
in figures 4, 5 and 7. In the absence of any spread in the transverse position and the angle of 
the beam, the charge spectrum would have been a step function of the applied plate voltage 
(V). It is the distribution of transverse position and the angular spread of the beam with 
respect to the exit slit that determines the shape of the spectrum. In the PSD mode of 
operation, the spectrum would be the exact map of the distribution of transverse position of 
the particles. In the voltage scanning mode, however, at each voltage point, a position 
integral, between the edges of the exit slit, of the position distribution of the beam on the exit 
slit plane is obtained. It should be noted that there is always some overlap in the integrals of 
the position distribution at successive voltage points. One needs, therefore, to deconvolute 
the voltage mode spectrum to find the fraction of particles of different charge states. 

Let us assume that there are particles of different charge states q's passing through the 
charge selector. Let N q be the number of particles of charge state q entering the selector at 
any plate voltage. In the voltage mode charge spectrum, F(V) is the total charge passing 
through the exit slit at the applied plate voltage V and collected by the Faraday cup 
behind the exit slit. 

If f q (V) be defined as the contribution of particles of charge state q to the total charge 
F(V) at the voltage point V, then 



) ainc) be the probability density function describing the charge independent 
distribution of the divergence parameters y[ nc and ai nc . 
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where 

W(s) = 9(s-si}- 6(s - s 2 ) (5) 

is the window function representing the exit slit with edges at si and 52, s 2 > s\. Here 

6(x) = 1 if x > 
= otherwise. 

Therefore, by change of variable u = C 2 (qV/E] 
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And 

ff q (V}dV = ^ /, (8) 
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where 
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The value of the integral / is independent of the charge state and depends only on 
^(yincjQiinc) and the geometrical constants of the charge selector. 

By similar change of variable, it is easy to show that (for a given energy) 

kN q , (10) 

g(qV) any function of qV : (11) 

k = constant indeoendent of a. f 12) 



suggested. The charge state fractions i.e. N q /Y^ q Nq can therefore be obtained directly 
from the areas under the peaks for different charge states. It may be noted that in PSD 
mode, the area under individual charge profile obtained by integrating the profile over the 
variable s would also be proportional to N q . 
In the special case of Gaussian distribution 
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it can be shown by explicit evaluation of the integral that 
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It is worth investigating two limits of the above results. 




I. Narrow exit window: When (s 2 s\ ) <C . o^, 4- <j^ 
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That is, f q (V] is a Gaussian with peak at 
s\ +s 2 E 



and 



FWHM = 2.35 



(14) 
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(16) 
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It may be noted that the peak height of f q (V) is proportional to qN q with the constant of 
proportionality being equal for all the charge states. So the peak heights can also be used 
to find the charge state fractions. 



II. Wide exit window: When (s 2 - s\) > 

f q (V) = qN q at V around V q = 
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5.1 Experimental setup 

The performance of the charge selector was checked using a setup discussed in our earlier 
papers [1,2]. Briefly, a highly collimated beam of energetic heavy ions in a given charge 
state obtained from the 14MV pelletron accelerator at Mumbai was made to fall on a 
target mounted inside an electrically isolated scattering chamber. In some cases the ions 
were further stripped using a post accelerator foil stripper [3]. The targets used were thin 
carbon foils of varying thicknesses mounted on a target wheel [4]. The scattered beam at 
entered the charge selector through an isolated vertical 1.5 x 12mm 2 slit, mounted just 
before the plates. The charge state distribution of the scattered beam was measured by 
deflecting it as a function of the voltage on the plates. The deflected beam, transmitted 
through a horizontal slit 10 x 2mm, was collected on a Faraday cup. A beam profile 
monitor, placed before the scattering chamber, was used to normalize the charge of the 
incident beam by measuring the number of emitted electrons from the BPM cathode wire. 
All the charge collected were digitized using AD converters which were specially desig- 
ned and made such that they could be used even when the incident beam intensity is very 
low in the measurement. Since the AD converters were sensitive to only one polarity, a dc 
offset was helpful in averaging out small negative pickup by the electrometer. The charge 
spectra as a function of the plate voltage were recorded several times for a particular 
incident charge state and energy. They were normalized for a fixed incident charge for 
further analysis. 
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by direct acceleration or by allowing the beam to go through a post-stripper [4] and 
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Figure 5. Normalized charge spectrum as in figure 4 with 96MeVF 8+ ions incident 
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Figure 6. Charge state fractions obtained from the spectra taken with the charge 
state selector for F 8+ ion incident on 30 ng/cm 2 C foil as a function of energy of the 
projectile. The solid curves are results reported by Shima et al [5]. 



subsequently selecting the desired charge state by a switching magnet. For obtaining 
several charge states to go through the selector the beam was passed through a set of 
carbon foils of varying thickness between 1.5 and 30 ug/cm 2 . Any one of the carbon foil 
could be brought onto the path of the beam by remote control. Some positions on the 
holder were left blank in order to allow ions in a definite charge state to enter the selector. 
Typical charge profiles obtained at the Farady cup with a 5 + or 8 + charge state of the 
incident beam are shown in figure 4 for 60 MeV F ions. The latter was obtained with the 
help of the post stripper foils. The observed line shape is Gaussian. The resolution in VS 
mode, AV/V where AV is the FWHM of the profile and V is the voltage at which the 
peak appears, was found to be ~ 5.3% in the present geometry for both the cases 
suggesting that the inherent beam quality (emittance, focussing and stripping elements) to 
be the deciding factor. In cases where the beam enters the selector through a carbon foil 
placed before the selector, because of the Coulomb interaction, the beam is split into 
several charge states. The charge state fractions thus obtained depend not only on the 
incident energy and the atomic number of the ion beam but to a significant extent on the 
thickness of the carbon foil used until charge equilibrium is attained (before equilibrium 
is reached these fractions may also depend on the incident charge state). This is nicely 
illustrated in figure 5 which shows the observed charge spectra for several thicknesses of 
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Figure 7. Normalized charge spectra, as in figure 4, with a beam of: (a) 24 MeV F 8+ 
(obtained using gas stripper) incident on 3 (ig/cm 2 C foil, (b) 56 MeV F 8+ (obtained 
using foil stripper) incident on 3|j.g/cm 2 C foil and (c) 56MeVF 8+ (obtained using 
foil stripper) incident on 30 ug/cm 2 C foil. 



cached. From a fit to the spectra with Gaussian line shape the charge state fractions were 
letermined (see IV). We have shown in figure 6, as a function of the incident beam 



with that of Shima et al. However, as seen from the higher^energy aaia. i is 
that charge equilibrium is not achieved with the 3()|ag/enr C foil used in ". in 
presently observed deviations at higher energies are expected to disappear had ue up 
thicker foils. Such measurements would result in an increase of the i r traction with .i 
corresponding decrease in the 8+ and 7+ fractions. The effect of the increased dt\eiveiKv 
of the beam due to the thickness of the C foil arising mainly due to multiple sc.ttii'iHr> 
can be seen in figure 7 for F ions of 56MeV. A comparison of the charge spectra taken 
with 3 and 30u.g/cm 2 foils (figures 7b and 7c) indicates the broadening t.I tin- !nu. 
caused mainly by Rutherford scattering. This is not so obvious at higher eneii'\ .i 
96MeV (figure 5) where this contribution is small. Figure 7a shows the eluufr sjuviiuiii 
with 24MeVF 8+ obtained using gas stripper in the terminal of the acceleratoi to !u\- a 
beam with smaller divergence. The charge spectrum taken with a 3 nn/enr' f Iml IN -.nil 
very well resolved even though one expects larger line-widths due to Kuiheitm.l 
scattering because of the still lower energy suggesting that the beam quality brfnjr rhr 
scattering chamber is the principal factor deciding the resolution. 

6. Conclusion 

An electrostatic charge selector for charge slate analysis has turn eonstnu'teil ,ui<! a 
performance tested. The data obtained using this selector on charge stale tt.u lions . } 
ions after interacting with C foil are found to he in good agreement with ihr puNi'-in*.! 
data of Shima et al [5]. The effect of beam divergence on the experimental sjuvUal lint- 
width is well parametrized by the analytical expressions presented ho re. 
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Abstract. Absolute total electron scattering cross sections for carbon dioxide have been measured 
at low electron energies using a photoelectron source. The measurements have been carried out at 
27 electron energies varying from 0.9 1-9. 14 eV with an accuracy of 3%. The cross sections 
obtained in the present experiment have been compared with other measurements and theoretical 
computations. 
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1. Introduction 

In recent years, there has been considerable interest in the measurement of total electron 
scattering cross sections of carbon dioxide, especially at low electron energies. Also, a lot 
of activity in the theoretical study of e-CO2 collisions has started in the recent past. One 
of the reasons for such an interest lies in the fact that CO2 plays an important role in the 
study of planetary atmospheres, laser fusion plasma and gas lasers. Also, this is one 
polyatomic molecule for which accurate low energy electron scattering cross section 
calculations are available. 

The electron scattering cross section measurements for carbon dioxide have been 
carried out using various techniques. This includes the work reported by Briiche [1], 
Ramsauer and Kollath [2], Szmytkowski and Zubek [3], Ferch et al [4], Hoffman et al 
[5], Sueoka and Mori [6], Szmytkowski et al [7] and Buckman et al [8]. These 
measurements show an increase in cross section with decreasing electron energy below 
1.8 eV and at around 70 meV, the cross sections have been found to exceed the value at 
the maximum of the 2 II H resonance at around 3.8 eV by a large factor [4]. The cross 
sections have also been theoretically computed by Morrison et al [9] using coupled 
channels procedure together with the adiabatic nuclei approximation. Comparison of the 
measured and the theoretically computed cross section values shows a general agreement 



helium, neon [10], argon, krypton, xenon [11], molecular hydrogen [12] and molecular 
oxygen [ 13] at electron energies from 0-10 eV. In this paper, we present the cross sections 
for carbon dioxide in the electron energy range between 0.91 and 9.14eV at an electron 
energy resolution of 45 meV. 



2. Experimental set up 

The experimental set up used for the study, the method for analysis of the data and error 
analysis have been described in detail previously by Kumar et al [10] and Subramanian 
and Kumar [11]. In brief, the experiment consists of measuring the intensities of the 
peaks in the photoelectron spectra of the source gas such as argon, krypton, xenon or 
carbon dioxide itself. The photoelectrons thus produced by the source gas are scattered by 
the target gas, the cross section of which is to be determined. Each peak in the photo- 
electron spectrum of the source gas provides one point in the electron energy scale. 

Electrons are produced by photoionization of the source gas in a small ionization region 
located close to the electron energy analyser. Vacuum ultraviolet photons from strong 
atomic emission lines are allowed to impinge on the source gas which is photoionized and 
electrons of different energies corresponding to the states of the ions are produced. With 
different combination of photons of three different wavelengths (He 1 58.4, NeI73.6 and 
74.4 nm) and argon, krypton, xenon and carbon dioxide as source gases, a large number of 
electron energies are accessible from 0.91 to 9.14eV. The resonant emission lines He I 
(58.4 nm) and Nel (73.6 and 74.4 nm) are produced by microwave discharge of helium and 
neon respectively. The intensity of the photon beam is monitored using a beam splitter 
which consists of a high transparency wire mesh mounted at an angle of 45 to the beam 
axis. More than 90% of the light is transmitted and the rest is reflected towards a light 
detector which includes a perspex light pipe coated at the front end with a scintillator, p- 
terphenyl, and an EMI 6199 photomultiplier. Any change in the light intensity could thus 
be monitored using the beam splitter, and the change in intensity, if any, during the 
experiment could be incorporated in the results. 

The transmitted photon beam was allowed to pass through a large number of circular 
baffles before emerging into a small region where photoionization of the source gas takes 
place. The pressure of the source gas was kept low to avoid any electron scattering by the 
source-gas atoms/molecules. The photoelectrons thus produced were energy analysed 
using a cylindrical mirror analyser (CMA) and were detected by a channeltron (4029, 
Galileo Electro-Optics Corp., USA) operated in the counting mode. After proper ampli- 
fication, the signals were stored in a microprocessor-controlled multichannel analyser 
operated in the multiscaling mode. 

The CMA has a slit-to-slit focusing geometry and the diameter of the inner and outer 
cylinders were 10.2 and 25.4cm respectively. The photoelectrons produced in the ioniza- 
tion region were allowed to enter the first slit of the analyser at an angle of 5444' 3. 



iu ueiuie me jj angie over trie ZTT geometry, was z.o / cm. wimoui any targei gas, 
the amplitudes of the photoelectron peaks produced by the source gas were monitored at a 
fixed but low source gas pressure in the ionization region. The target gas was then 
introduced and if the source gas and target gas were not the same, the amplitudes of the 
photoelectron peaks would decrease with increasing target gas pressure. This decrease 
was attributed to electron scattering by target gas species whose number density increases 
with increasing pressure. When the source and target gas were the same (e.g. in case of 
carbon dioxide), the amplitudes of the photoelectron peaks would first increase with 
increase of source/target gas pressure and after a certain gas pressure was attained, the 
amplitudes would start decreasing with subsequent increase in gas pressure. To begin 
with, the increase in peak amplitudes would be due to the increase in photoelectron 
production rate followed by small electron scattering, but the subsequent decrease in 
amplitudes was because of scattering alone. The absolute gas pressure in the ionization 
and scattering regions was measured using an MKS capacitance manometer (head 
310MH-1). The details about pressures in different regions of the instrument have been 
given in a previous publication [10]. 

3. Method and error analysis 

Using the technique described above, electron scattering cross sections for carbon dioxide 
have been measured at 27 energies. When source and target species are different, the 
electron scattering cross section a can be evaluated by using the following equation: 
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Here l e \ and l e i are the amplitudes of the photoelectron peak at two different target gas 
pressures P\ and P 2 , /AOI and /AOI are the incident photon intensities monitored by the 
beam monitor during the experiment at two pressures, l\ is the distance from the centre of 
the beam splitter to the circular aperture covering the ionization region, h is the distance 
from the same aperture to the actual ionization region defined geometrically by different 
slits in the accelerator region and cylindrical mirror analyser, x is the scattering path 
length, a is the ratio of the pressures of the target gas outside and inside the ionization 
region and k is the photoabsorption coefficient of the target gas at a particular photon 
wavelength (figure 1). The total photoabsorption coefficients for carbon dioxide at the 
three incident photon wavelengths have been taken from Samson and Yin [14]. The ratios 
of /AOI and /A02 were determined from the beam splitter as described previously [10]. 

When source and target gas species are the same, the electron scattering cross sections 
could be evaluated by the equation given below: 
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Figure 1. Schematic diagram showing ionization as well as the electron scattering 
region in the cylindrical mirror analyser of the photoelectron spectrometer. 



present experiment was estimated to be 3%. The electron scattering cross sections for 
carbon dioxide were measured using equations (1) and (2) respectively in the two cases 
described above. In both these equations, all parameters I e \ , I e i, /AOI i Ao2 P\ and PI could 
be determined experimentally and cross sections could be calculated. 

The problem of inadequate discrimination against forward scattered electrons is a 
difficult one. Such an inadequacy is due to the inability of the experimental system to 
discriminate low angle scattered electrons which, in general, have large amplitudes. 
Attempts to estimate the magnitude of this effect on the measured cross sections in a 
certain experimental set up becomes much more difficult as almost no experimental 
values of elastic and inelastic differential scattering cross sections at very low angles are 
available. An attempt, however, has been made in this direction for our experimental set 
up using a numerical simulation method. 

In the present experimental set up, the geometry of the apparatus allows a full cone 
angle of 3 beam of electrons into the analysing region (figure 1). Taking into account the 
forward scattering of the electrons near the exit slit of the scattering region, it would be 
safe to assume that the overall angle of emergence is increased by a factor of two, giving 
the full cone angle of 6 around the nominal angle of projection of 54 '44'. 

In order to study the problem of inadequate discrimination against forward scattered 
electrons, a Monte Carlo method was developed. The electron trajectories originating 
from the ionization region and terminating at the detector were numerically studied using 
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are produced, (ii) random angle of projection within the allowed cone of ejection 
determined by the geometry, (iii) spread in energy of the photoelectron determined by the 
thermal broadening, (iv) place nearby the first slit A (figure 1) where a chosen fraction of 
the emitted photoelectrons are undergoing forward scattering and (v) deflection in angle 
within a specified cone due to the forward scattering. Assuming cylindrical symmetry of 
the spectrometer, the analysis was carried out in the in-plane geometry. In the absence of 
values of differential cross section available at small angles, the forward scattering cross 
section at the exit slit was arbitrarily chosen as 100 A 2 per steradian, which is about 6 
times larger than the peak value of the measured total cross section in the present work. 
The discrimination due to geometry was computed using normal trigonometrical 
relations. The discrimination due to the range differences are computed with the standard 
equation for range of electron trajectory in the analysing field of a slit-to-slit focussing 
type CMA given by 

L~r\ x \/47rfccos 0exp( sin 2 0)erf (Vfcsin 0). (3) 

Equation (3) has been derived by Karras et al [15]. The error function erf (^/k sin 9} is 
defined as 
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where x = \/ksin6. The analyser constant, k = E/eVlnfa/ri], E is the energy of the 
electron in eV, e is the electronic charge, V is bias potential applied to the analysing field 
electrode, is the angle of projection and r\ and r 2 are the radii of the inner and outer 
cylinders of CMA. 
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standard deviation values. Finally, similar data were generated for different angles of 
emergence for forward scattering giving the full cone angles of 2,4,6, 8 and 10. 

The result shows that the average value of error in cross section due to lack of 
discrimination is more or less independent of Ae/E in the chosen range of 0.2 to 0.6%. 
However, it is very much dependent on the overall cone angle of emergence. The results 
of the present study are shown in figure 2, where the error in cross section value is plotted 
against the overall cone angle of emergence. 

For our experimental system where the overall cone angle of emergence has been 
estimated to be 6, the maximum error due to inadequate discrimination against forward 
scattering has been found to be about 0.5%. 



4. Results and discussion 

The total electron scattering cross sections for carbon dioxide as measured in the present 
experiment are shown in figure 3 for electron energies ranging from 0.91 to 9.14eV 
along with error bars at four energies only. Also shown in the figure are the measured 
cross sections reported by Ferch et al [4], Hoffman et al [5], Buckman et al [8], and 
Szmytkowski et al [7] and theoretically computed cross sections by Morrison et al [9]. 
Transmission measurements made with an electron time-of-flight spectrometer were 
reported by Ferch et al [4] in the low electron energy region (0.07 to 4.62 eV) with an 
energy resolution of lOmeVat 0.07 eVand 250 meVat 3.8 eV respectively. The measure- 
ments of Hoffman et al [5] were carried out using a beam transmission technique at 
electron energies from 2-50 eV at an energy resolution of 140 meV. The retarding 
potential time-of-flight technique was used by Sueoka and Mori [6] and a time-of-flight 
technique was also used by Buckman et al [8] to measure cross sections at electron 
energies from 5 to 403 eV and 0. 1 0-5.0 eV with an energy resolution of 300 and 250 meV 
respectively. A linear transmission method was employed by Szmytkowski et al [7] for 
electron energies varying from 0.5 to 80 eV at the energy resolution of 50 meV. The cross 
sections were also measured using the Ramsauer technique by Briiche [1] and Ramsauer 
and Kollath [2] and the normalized values of cross section were reported by Szmytkowski 
and Zubek [3] using a linear transmission method in the electron energy 1-49 eV, 0.2- 
2 e V and 2-8 eV respectively. The theoretically computed cross sections were reported by 
Morrison et al [9], Morrison and Lane [16], Collins and Morrison [17], and Lynch et al 
[18]. A theoretical coupled channels investigation was reported by Morrison et al [9] for 
incident electron energies from 0.07 to 10.0 eV. The fixed nuclei approximation was made 
with the molecule in the ground state and the nuclei frozen at their equilibrium positions. 
The static-model exchange-polarization interaction potential was used to calculate total 
integrated cross sections. The improvement in the method was carried out by Collins and 
Morrison [17] at electron energies ranging from 0.1 to 2eV by using orthogonalized 
model potential. Lynch et al [18] calculated the integrated elastic electron scattering cross 
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Figure 3. Total electron scattering cross sections for carbon dioxide as a function of 
incident energy from 0.91 to 9.14eV obtained by various researchers. The symbols 

denote: o, Buckman et al [8]; V, Hoffman el al [5], Morrison et al [9]; Q, 

Szmytkowski et al [7]; A, Ferch et al [4]; , present work; , cross 

section curve after proper curve fitting (present work). 



section values reported by Briiche [1], Ramsauer and Kollath [2], Szmytkowski and 
Zubek [3], Morrison and Lane [16], Collins and Morrison [17] and Lynch et al [18] have 
not been included in figure 3. 

For the sake of discussion, the cross section curve (figure 3) in the total electron energy 
range from 0-9 eV has been divided into three regions. These include energy regions 
below 2 eV, between 2 to 5 eV and between 5 to 9 eV. In the first energy region, the cross 
section has been found to decrease rapidly with increasing energy followed by a sharp 
increase at higher energies. This rapid decrease is mainly due to s-wave scattering and 
thus can be explained by a nearly bound / = state in the e-COi potential [9, 19]. The 
existence of such a virtual state has been verified both experimentally by Kochem et al 
[20] as well as theoretically by Estrada and Domcke [21]. At energies below 1.1 eV, there 
is an excellent agreement in cross section values reported by Buckman et al [8], 
Szmytkowski et al [7], Ferch et al [4], Morrison et al [9] and those obtained in the present 
experiment. But the agreement worsens at increasing energies. There is a large 
discrepancy between the results shown in figure 3 compared to values reported by Lynch 



work and by Szmytkowski et al [7] to be 6. 13 and 5.84 A 2 at 1.88 and 1.90 eV respectively 
whereas Ferch et al [4], Hoffman et al [5], Buckman et al [8] and Morrison et al [9] have 
reported the cross section values of 5.85, 6.02, 5.84 and 6.21 A 2 at the minima of around 
1.8, 2.0, 2.0 and 2.0 eV respectively. It may be pointed out here that the minima stated 
above are not those directly measured by the respective groups, but rather are the results 
obtained by fitting curves to the values measured by them. Also, it may not be possible to 
have an exact energy minima in case of Buckman et al [8] and Hoffman et al [5]. In the 
first case [8], the energy steps used are large (0.5 eV) and therefore, elude finding the 
exact minimum; in the second case [5], there is no data available below 2.0 eV. A 
comparison of all these values shows that the largest disagreement in cross sections at the 
minimum is about 6.3%. 

The results obtained in the present experiment in the energy region 2-5 eV (figure 3) 
shows a sharp increase in cross section up to the well known maximum centered around 
3.8eV followed by a sharp decrease up to 5eV energy. The peak around 3.8eV was 
attributed to the formation of a shape resonance [22]. The temporary state involved was 
the 2 II U of the CO^T ion having a short lifetime. In general, the shape of the resonance 
peak has been faithfully reproduced by all the work shown in figure 3 but the cross 
section values at the maximum of the resonance peak and also at the wings around the 
peak show some discrepancies. To have a clear comparison, the cross section curves 
reported by Ferch et al [4], Buckman et al [8], Hoffman et al [5], Szmytkowski et al [7], 
and Morrison et al [9] along with the curve obtained in the present experiment have been 
made to undergo proper curve fitting procedure. The cross section curve after proper 
curve fitting to the results obtained in the present experiment has also been shown in 
figure 3. This has been carried out using an appropriate mathematical function in the 
library of large number of such functions available in the Tablecurve-2D (Jandel 
Scientific) software. After proper fitting, the values of electron energy and cross section 
at the peak of the shape resonance, width of the peak at 75% of the maximum value and 
the slopes of the peak at 75% value on both wings of the peak have been obtained and 
these are given in table 1. The values of 75% was chosen as in most of the curves, it was 
difficult to get full width at half maximum. Also given in the table are the values of 
energy resolution used by researchers in their respective experiments. From the table, it is 
clear that there is a good agreement in the energy values obtained at the peak of the 
resonance but the cross section values show a large discrepancy between the present 
results and other investigations. There is large overestimation of cross section of 30.88 A 2 
by Morrison et al [9] and this is most likely due to the neglect of nuclear motion in the 
calculation [8]. The cross section values given by Ferch et al [4] and Buckman et al [8] 
are 14.2 and 5% lower as compared to the cross section reported in the present work 
whereas the cross section value (not given in table 1) reported by Sueoka and Mori [6] is 
almost 30% smaller. There is a fairly good agreement in the cross section value reported 
by Szmytkowski et al [7] and Hoffman et al [5] as compared to that given in the present 
experiment. It is expected that higher the energy resolution used in the experiment, higher 



given in table is the energy resolution used in different experiments. 
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other within the experimental uncertainties except possibly those of Ferch et al [4]. It 
appears that such a behaviour cannot be explained by the above argument alone but it is 
felt that the above behaviour could be explained by the combined effect of narrower 
electron energy width and better discrimination against small angle scattering in different 
experiments. In the present experiment, it has been shown that error in cross section 
values due to lack of discrimination against low angle scattering is about 0.5%. Hoffman 
et al [5] have pointed out that their measurements may be 1% too low at 4eV and less 
than 1% too low at 10 eV because of inadequate discrimination against forward scattered 
electrons. In a subsequent paper, Buckman and Lohmann [23] have discussed about 
the problem of low angle scattering in their experimental system and estimated this 
error to be less than 0.1% at about 2 eV electron energy. The small angle contribution 
was estimated to be below 0.5% at all electron energies of interest in the experimental 
system reported by Szmytkowski et al [7] where as no such estimate was reported by 
Ferch et al [4]. 

Also given in table 1 are the AE values at 75% of the peak value of the shape resonance 
and modulus of the slopes at the two sides of the peak at 75% peak value. The values of 
the modulus of the slopes on the rising and falling sides of the peak as given by different 
experimental groups show that the resonance peak is more or less symmetrical around the 
maximum in the present work as well as in the experiment reported by Szmytkowski et al 
[7] whereas the experiments by Ferch et al [4], Buckman et al [8] and Hoffman et al [5] 
show that the resonance peak is asymmetrical in nature. This may be attributed to the 
energy resolution used in the experiment. It has been found (table 1) that higher the 
energy resolution, more symmetrical is the resonance peak around the maximum. 

In the energy region from 5-9 eV, cross section values reported by the three experi- 
mental groups only ([5,7] and present work) have been shown in figure 3 along with the 
theoretically computed values of Morrison et al [9]. Unfortunately, Ferch et al [4] and 
Buckman et al [8] have not carried out these measurements beyond 4.62 and 5.0eV 
respectively. The general trend in this energy region as suggested by different research 



Wavelength 
A (A 2 ) 


Source Photo-ion 
gas state 


Electron 
energy 
(eV) 


Cross sect 
(A 2 ) 


744 


Ar 2 P 3/2 


0.91 


9.38 


736 


Ar 2 P 3/2 


1.09 


7.96 


584 


CO 2 C 2 S+ 


1.80 


6.01 


744 


Vr 2p 
JtV.1 r 1 /2 


2.00 


6.19 


736 


Kr 2 P 1/2 


2.18 


6.84 


744 


Kr 2 P 3/2 


2.66 


7.83 


736 


Kr 2 P 3/2 


2.85 


7.84 


744 


co 2 -^ 2 n,, 


2.90 


7.44 


584 


/T~\ D2v+^i, t\ 

1 \-J 2 LJ ' ' (Is J ) 


2.96 


8.15 


736 


co 2 x 2 n 


3.07 


8.80 


584 


r*r\ o2v i +/',i c\\ 
L-v_/ 2 D 2-i ' \y = U ) 


3.12 


9.74 


744 


Yo 2p 

Ae 'i/2 


3.23 


10.65 


584 


CO 2 A 2 H M (zx = 4) 


3.31 


10.62 


736 


Xe 2 P 1/2 


3.41 


12.89 


584 


CO 2 A 2 n, ( (y = 3) 


3.46 


13.74 


584 


CO 2 A 2 II M (zx = 2) 


3.61 


15.75 


584 


co 2 A 2 n i( (i/ = i) 


3.75 


15.68 


584 


C0 2 A 2 Ti u (v = 0) 


3.89 


16.22 


744 


Xe 2 P 3/2 


4.59 


8.95 


736 


Ae * 3/2 


4.77 


9.21 


584 


Ar 2p 

f" '1/2 


5.28 


6.71 


584 


Ar 2p 
z 3 / 2 


5.46 


7.18 


584 


Kr 2 P 1/2 


6.55 


7.73 


584 


&2p 
'3/2 


7.22 


8.02 


584 


PO^ Y^TT 
V *-' 2 A iig 


7.44 


8.85 


584 


Y 2p 

''1/2 


7.78 


8.88 


584 


Yo 2p 

Ac ' 3 /2 


9.14 


10.43 



et al [7] and Hoffman et al [5] occurs almost at the same energy of about 5.4eV 
different cross section values of 6.9, 7.9 and 8.5 A 2 respectively. Morrison et a 
reported a minimum value of 7.1 A 2 at the electron energy of 6.0eV but the increa 
cross section values was rather slow at higher electron energies. At electron ene 
between 5.5 and 9eV, rate of change of cross section per unit electron energy has 
found to be almost similar in the three cases, but absolute values of cross section 
quite different. At the two typical energies 7.5 and 9.0 eV, the cross section v 
reported in the present work are lower by 10% and 7% as compared to those give 



the corresponding electron energy. 
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Abstract. We have tried to understand the field dependence of magnetization of high temperature 
superconductors in the light of phenomenological theory. Especially, the field dependence of 
dM/dln5 of polycrystalline Bi(2212) is understood by incorporating the overlap of vortices in the 
London theory. 
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The magnetic properties of high temperature superconductors (HTSCs) in magnetic fields 
are very unusual near the transition temperature T c and are different for different HTSCs. 
Near T c , the magnetization (M) of Y(123) is proportional to (T - T*) 2 , where T* is 
almost field independent and close to T c [1,2]. It has also been observed that the 
magnetization of Y(123) is approximately logarithmic in B a in the reversible region and 
can be understood by using 3D description of the London theory [3-5]. It has also been 
applied to understand the temperature dependence of specific heat in magnetic field [6], 
whereas, the field dependence is understood by incorporating the overlap of vortices in 
the London theory [7]. Incorporating the overlapping of vortices in the London theory the 
magnetization for 3D superconductors can be written as [7] 

B-H fo Bg\ fof* , B C2 



This is similar to the result obtained from the variational model [8] in the London limit. The 
temperature dependence of penetration depth and the upper critical field is given by [9] 

A(T)= A and B c2 = B c20 (l - f) 2 ", (2) 

where t = T/T C . The v is 1/2 in the mean field region and 2/3 in the critical region. 
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Figure 1. The magnetization (in units of (f)Q/32'n 2 X 2 ) vs b. 



which are in between that obtained from the usual London theory and the variational 
approach proposed by Hao et al [7] and the magnetization is approximately logarithmic 
in B a . Near B c i or close to T c , the logarithmic term in eq. (1) can be expanded as 
(1 - (B a /B c2 ). So the magnetization behaves as (B c2 - B a ) 2 or (T c - T) 2 . If the upper 
critical field is large, the change on T c with magnetic field is little as observed 
experimentally [1,2]. 

The magnetic properties of Bi(2212) are different from Y(123). In the intermediate field 
range, the 4?rM vs \nB a curve of Bi(2212) single crystal is not linear as predicted from 
the London theory. The variation in the case of single crystal [10, 1 1] is different from the 
polycrystalline case [12]. According to the mean field theory developed by Koshelev [13] 
(from eq. Cll), the field dependence of the slope dM/dlnB can be written as 



dM 



327T 2 A 2 (r) 



Ba 



with 77 = 0.35. 



Whereas, according to the vortex overlapping mechanism [7], it is given by 



(3) 



(4) 



Shown in figure 2 is the dM/dlnZ? versus B a . The dashed line represents the Koshelev's 
result and the solid line represents the vortex overlapping result. It is noted from this 
figure that both mean field theory and vortex overlapping mechanism predicts different 




Figure 2. dW/dlnB. (in units of 

taken B c2 = 10 and for dashed line 5 c2 = 20 T. 



For the solid line we have 
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fit according to eq. (4) with parameters A = 



and B a = 



Genoud et al [\2] have claimed that the polycrystalline data can be fitted to 

100 atV= 55 K T t s e values are reasonaHy in good agreement 
*at the field dependence of dM/dln* can be 
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